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Abstract. Spin Density Matrix Elements (SDMEs) describing the angular distribution of exclusive 
p° electroproduction and decay are determined in the HERMES experiment with 27.6 GeV beam 
energy and unpolarized hydrogen and deuterium targets. Eight (fifteen) SDMEs that are related 
(unrelated) to the longitudinal polarization of the beam are extracted in the kinematic region 
1 GeV^ < Q2 < 7 GeV^ 3.0 GeV < < 6.3 GeV, and -t < 0.4 GeV^ Within the given 
experimental uncertainties, a hierarchy of relative sizes of helicity amplitudes is observed. Kinematic 
dependences of all SDMEs on and t are presented, as well as the longitudinal-to-transverse 
electroproduction cross section ratio as a function of Q^. A small but statistically significant 
deviation from the hypothesis of s-channel helicity conservation is observed. An indication is seen of 
a contribution of unnatural-parity-exchange amplitudes; these amplitudes are naturally generated 
with a quark-exchange mechanism. 
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1 Introduction 

In exclusive production of vector mesons such as p, uj 
or ^ from deep- inelastic lepton scattering (see Fig. [T|) , 
measurements of angular and momentum distributions 
of the scattered lepton and vector meson decay prod- 
ucts allow one to study the production mechanism and, 
in a model-dependent way. the structure of the nu- 
cleon. 

For more than 40 years, many basic features of vec- 
tor meson production by a virtual photon have been 
successfully explained in terms of the Vector Meson 
Dominance (VMD) model [HE] - In this model, the vir- 
tual photon fluctuates into a vector meson whose inter- 
action with the nucleon could be described, for exam- 
ple, using Regge phenomenology. More recently, in the 
context of perturbative Quantum Chromo-Dynamics 
(pQCD), exclusive meson production at sufficiently 
large values of the photon virtuality and the invari- 
ant mass of the photon-nucleon system W is assumed 
to be dominated by so-called handbag-diagrams (see 
Fig. [5]) that involve various non-pcrturbative nucleon 
structure functions, known as Generalized Parton Dis- 
tributions (GPDs) [31I41I51I6]. 

In pQCD, the common model of the production 
of vector mesons at high and W can be consid- 
ered as three consecutive steps [7]: i) dissociation of 
the virtual photon into a quark-antiquark (55) pair, 
ii) scattering of the pair on a nucleon (nucleus), iii) 




N(p) t N'(p') 

Fig. 1. A generic f-channel exchange process for 7*7V — > 
p^N' . Each particle's four-momentum is denoted in paren- 
theses. 

formation of the observed vector meson from the qq- 
pair. (A full quantum mechanical treatment includes 
all possible time orderings, which may be more impor- 
tant at lower energies.) The interaction of the gg-pair 
with the nucleon can proceed via two distinct mecha- 
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nisms. The first one, two-gluon exchange, is described 
by the Fcynman diagram shown in Fig. [2^. This pro- 
cess transfers the same quantum numbers as pomeron 
exchange in the Regge picture, and is anticipated to 
exhibit a similar phenomenology. The second mecha- 
nism is described by the exchange of a gq-pair, also 
possibly with additional gluons connecting them, and 
is called quark-exchange (Fig.[2lD). The corresponding 
process in Regge phenomenology [5] is the exchange 
of "secondary" reggeons, such as p, w, /2 and 02 in 
the case of natural-parity exchange (NPE), in which 
the spin J and parity P associated with the reggeon 
trajectory are = 0+, 1~, 2+, or tt, ai, hi mesons 
with = 0~, 1"'', ... in the case of "unnatural-parity" 
exchange (UPE). In the GPD formahsm, NPE (UPE) 
processes are described by H and E [H and E) GPDs. 
In the intermediate energy range of the HERMES ex- 
periment (3 GeV < W < 6 GeV) and the moderate 
values of photon virtuaHty (1 GeV^ < < 7 GcV^) 
both Regge phenomenology and pQCD may be applied 
to describe exclusive vector meson production. The in- 
terpretations they offer of the experimental data are 
often complementary, although not necessarily consis- 
tent. 

The main focus of this work is on the measurement 
of Spin Density Matrix Elements (SDMEs) of the p° 
meson, which describe the distribution of final spin 
states of this produced vector meson. These elements 
depend on amplitudes for the angle- and momentum- 
dependent transition processes between initial spin 
states of the virtual photon and final spin states of the 
produced vector meson. The values of SDMEs serve 
to establish the hierarchy of helicity amplitudes that 
are commonly used to describe exclusive p'^ produc- 
tion. In this way the relative importance of the vari- 
ous 7* transitions is revealed. Two main order- 
ing principles are observed in vector meson Icptopro- 
duction, s-channel helicity conservation (SCHC) and 
the dominance of NPE over UPE mechanisms. SCHC 
implies that only 7* p° transitions with the same 
helicities of virtual photon and p'^ occur in the reac- 
tion when considered in the "hadronic" center-of-mass 
frame (defined below). These concepts apply both in 
the reggeon-exchange picture and in pQCD. In par- 
ticular, we note that a signal of UPE is evidence of 
quark-antiquark exchange (Fig. Wp)^ as the pomeron 
has natural parity. 

At high energies pomeron exchange dominates, and 
secondary-reggeon exchanges with natural parity are 
suppressed by a factor ~ M/W [8] in their amplitudes; 
AI is an energy scale in Regge phenomenology cho- 
sen to be equal to the nucleon mass. Also suppressed, 
by a factor ~ {M/WY [H], are the most important 
unnatural-parity exchanges mediated by tt, ai, and bi 
reggeons. Therefore substantial UPE contributions can 
be expected only at lower values of W . 

In the pQCD framework, the leading-twist contri- 
bution describes the transition of longitudinal photons 




Fig. 2. Examples of a) a two-gluon exchange diagram and 
b) a quark-exchange diagram, shown for the lowest order 
in the strong coupling constant as ■ 



to longitudinal vector mesons, which is s-channel he- 
licity conserving and corresponds to natural-parity ex- 
change. As it is not agreed how strongly the various 
other contributions are suppressed at a given energy, 
measurements of SDMEs in the HERMES kinemat- 
ics help to distinguish these contributions and are of 
particular interest. Non-conservation of s-channel he- 
licity in exclusive p'^ production was already observed 
at collider energies pi fTnUTT] . At lower energies it was 
observed at HERMES [12], and for exclusive oj pro- 
duction at CLAS [13] . 

At sufficiently large values of W ^ experiments are 
typically sensitive to partons that carry small nucleon 
momentum fraction x, where the parton density in the 
nucleon is dominated by gluons. High-energy data of 
HI and ZEUS [9]IT0]ITT][T4] are well described by two- 
gluon exchange. At lower values of W , larger values 
of X are probed, where the parton density in the nu- 
cleon receives significant contributions from quarks. 
Indeed, a contribution from the quark-exchange mech- 
anism has been suggested to be necessary to describe 
exclusive p° production at intermediate virtual-photon 
energies, as in the case of the HERMES data [TSlfTBl 
[T7lfT§] and corresponding calculations [l9 l l20 ] l2T ] [22] . 

In leptoproduction, the spin transfer from the vir- 
tual photon to the vector meson is commonly described 
by helicity amplitudes, from which SDMEs can be con- 
structed. The detection of the scattered lepton and the 
vector meson decay products allows one to reconstruct 
the full reaction kinematics and the three-dimensional 
angular distribution of the production and decay of 
the p^ meson. For an unpolarized or helicity-balanced 
lepton beam, the expression for this distribution con- 
tains a set of "unpolarized" SDMEs as coefficients. An 
additional set of "polarized" SDMEs, which appear in 
products with the beam polarization in the expression 
for the angular distribution with polarized beam, can 
be determined if information on the longitudinal polar- 
ization of the lepton beam is available [23l[24] . In a very 
recent new classification scheme of SDMEs [5S], also 
the cases of longitudinal and transverse target polar- 
izations are described. However, the analysis in this pa- 
per follows the representation introduced in Ref. 
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Early theoretical calculations |2] of SDMEs in vec- 
tor meson production were based on the VMD model. 
More recent calculations combining this model with 
pQCD models [7l[26ll27ll28l[29ll30] and with Regge phe- 
nomenology [3TII32] are mainly focused on the high- 
energy kinematics of the HERA collider data. A con- 
temporary account of the various theoretical approaches 
is given in Ref. [M]. Recent model calculations based 
on GPDs present SDMEs for both high and inter- 
mediate energies, considering first only two-gluon ex- 



In the laboratory (lab) frame, ?? is the scattering an- 
gle between the incoming and outgoing Icpton, whose 
incoming and outgoing energies are denoted by E and 
E' . The photon virtuality is given by: 



« AEE'srn^ -, 
2 



(2) 



which is positive in leptoproduction. In this equation 
the electron rest mass is neglected. The four-momenta 



change [33], and recently incorporating quark exchange [35 the target nucleon and of the recoihng baryon are 

denoted by p and p' , respectively, and both have rest 



In this analysis, the beam polarization is used for 
the first time in an SDME extraction, thereby making 
possible the determination of the additional 8 polar- 
ized SDMEs. The high-statistics data samples accu- 
mulated at HERMES in the years 1996-2005 on both 
hydrogen and deuterium targets are used to determine 

decay angle distributions with an accuracy superior 
to that of the previously published HERMES ^He data 
from 1995 [16] and of the preliminary HERMES results 
from hydrogen data collected in 1996-1997 [T2l[36] . 
The improved statistical accuracy permits the study 
of the nature of the exchange mechanism, and in par- 
ticular the testing of the hypothesis of s-channel he- 
licity conservation. The availability of both hydrogen 
and deuterium targets offers the possibility to search 
for significant contributions of secondary reggeon ex- 
change with isospin / = 1 and natural parity. 

The structure of this paper is as follows. The kine- 
matics, SDME formalism, and HERMES experiment 
are described in the next three sections. The analysis 
procedure including event selection and background 
subtraction is discussed in section 5. The extraction 
of the SDMEs from the data using a Monte Carlo 
based maximum likelihood method is described in sec- 
tion 6. The experimental results on SDMEs integrated 
over the entire observed kinematic region are pre- 
sented in section 7, and their kinematic dependences 
are shown in section 8. An indication of the contribu- 
tion of unnatural-parity-exchange amplitudes is dis- 
cussed in section 9. Contributions of helicity-flip and 
UPE amplitudes to the cross section are estimated in 
section 10. The ratio of longitudinal to transverse 
electroproduction cross-sections is presented in section 
11. The results are summarized in section 12. 



2 Kinematics 

Figure [1] identifies the kinematic variables of lepto- 
production, 



(1) 



where N{N') denotes the initial (scattered) nucleon. 
The four-momenta of the incoming and outgoing lep- 
ton are denoted by k and fc', the difference of which 
defines the four-momentum q of the virtual photon 7* . 



mass M of the nucleon. irrespective of target. 
The Bjorken scaling variable is defined a^ 



Xb 



2p-q 2Afi/' 



with 



p_q 

M 



lab 



E - E\ 



(3) 



(4) 



so that v represents the energy transfer from the in- 
coming lepton to the virtual photon in the labora- 
tory frame. The squared invariant mass of the photon- 
nucleon system is given by: 



= (q+pf ^M^ 



(5) 



The squared four-momentum transfer from virtual pho- 
ton to vector meson equals that between the momenta 
of the initial and final nuclcons or nuclei. 



t^{q-vf ^{p-p'f 



(6) 



where v is the four-momentum of the vector meson. 
The variables t, Iq, and 



t' ^t~tn 



(7) 



are always negative, where —to represents the smallest 
kinematically allowed value of —t at fixed v and Q^. 
In the photon-nucleon center-of-mass frame considered 
here, the condition t = to corresponds to the case 
where the momentum of the produced p° is coUinear 
with that of the 7*. Typically for exclusive processes 
at intermediate and high energies, \to\ is much smaller 
than \t\ and therefore t' w t. 

At very low t, the approximation —t' w holds, 
where vt is the transverse momentum of the vector 
meson with respect to the direction of the virtual pho- 
ton, i.e., the subtraction of to removes the contribu- 
tion of the longitudinal component of the momentum 
transfer. 



^ This kinematic observable is to be distinguished from 
the variable x of the quark parton model, which represents 
in the GPD formalism the average longitudinal momentum 
fraction of the probed parton in the initial and final states. 
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scattering-plane 




Fig. 3. Definition of angles in the process 7*A*' — > 
p°N' -K+n'N' [37]. Here <5 is the angle between the 
p° production plane and the lepton scattering plane in the 
"hadronic" center-of-mass system of virtual photon and 
target nucleon. and are polar and azimuthal angles of 
the decay vr^ in the vector meson rest frame. 



The variable e represents the ratio of fluxes of lon- 
gitudinal and transverse virtual photons: 



1-2/ 



lab I 

^ l + 2(l + ^)tan2| 

lab 



(8) 



with y = p ■ q/p ■ k = v/E. 

The "exclusivity" of p° production is characterized 
by the variable 



AE 



Ml 



lab 



2M 



Ex, 



(9) 



where Mx is the invariant mass of the recoiling sys- 
tem, Ey = v + t/{2M) is the energy of the exclusively 
produced p° meson, and {E^+ + E^- ) is the sum of the 
energies of the two pions. For exclusive vector meson 
production H]), Mx = M holds and hence AE = 0, 
given perfect detector and beam energy resolution. 
Angles used for the description of the process 
— > p°N' — > Tr'^Tr~N' are defined according to 
Ref. [37] and presented in Fig. [3l The helicity amp- 
litudes are defined in the "hadronic" center-of-mass 
system of virtual photon and target nucleon, where 
the Z-axis is directed along the virtual photon three- 
momentum q. The y-axis of the right-handed system 
is parallel to q x v. It is the normal to the p° produc- 
tion plane spanned by the three-momenta q and v, of 



the virtual photon and p'^-meson, respectively. The an- 
gle between the p^-production plane and the lepton- 
scattering plane in the "hadronic" center-of-mass sys- 
tem is specified by: 



cos^ = 



sin<? 



(q X v) • (k X k') 
|q X v| • |k X k'l ' 
[(q X v) X (k X k')] • q 
|qxv|.|kxk'|.|q| ■ 



(10) 



The angle (f) between the p^-production plane and 
p°-decay plane is defined by: 



cos ( 



sm = 



(q X v) ■ (v X p^+) 
|q X v| • |v X p^+l ' 
[(q X v) X v] ■ (p^+ X v) 
|(q X v) X v| • |p^+ X v| ' 



(11) 



where p^+ is the three-momentum of the positive de- 
cay pion in the "hadronic" center-of-mass system. 

The polar angle of the decay 7r+ in the vector 
meson rest frame, with the z-axis aligned opposite to 
the outgoing nucleon momentum P' and the y-axis 
parallel to Y and directed along P' x q, is defined by: 



cos = 



P' P, 



(12) 



where P7r+ is the three-momentum of the positive de- 
cay pion. 

Note that the relation between this notation and 
the notations of the so-called "Trento convention" [38] 
and Ref. |25] is: ^ = -^ps] = -0/i[38], 



6* = i9 



^[25], 



[25]- 



3 Formalism 

3.1 Helicity Amplitudes 

Exclusive vector meson leptoproduction ([l} is com- 
monly described by helicity amplitudes i^AvA'„;A.^Aiv , 
defined in the "hadronic" center-of-mass system of vir- 
tual photon and target nucleon [23] (see Fig. [3]) . He- 
licity indices and Ay describe the spin states of 
virtual photon and p meson, respectively, while Xx 
{X'x) is the helicity of the initial (scattered) nucleon. 
The helicity amplitude can be expressed as the scalar 
product of the matrix element of the electromagnetic 
current vector J*^ and the virtual-photon polarization 



vector e 



F: 



{vXv;p'X'xnXN)el^^\ (13) 



where e^^^^ describes the transverse and e^°^ the lon- 
gitudinal polarization of the virtual photon. The ket 
vector IpAat) corresponds to the incident nucleon and 
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the bra vector {v\v\p' describes the final state of 
the p'^ meson and scattered nucleon. The amplitudes 
depend on Q^, W and t. For convenience, these depen- 
dences may be omitted in the following. 
The amplitudes obey the relation [23] 

(_l)(Av-AW)-(A,-A«)^^^^^,^^^^^ (14) 

which is a consequence of parity conservation in the 
strong and electromagnetic interactions. 



3.2 Natural and Unnatural-Parity-Exchange 
Amplitudes 

A helicity amplitude F can be decomposed into an 
amplitude T for natural-parity exchange and an am- 
plitude U for unnatural-parity exchange: 

-^Av A^;A-^ Aat ^AvA^;A^Aiv ~^ ^AvA^;A^Ajv' (-^^) 



with 



Av A^ ; A^ Aiv 



^(^Av-A^;A^A„ + (^1) '^''"'''^^^-AvA;^;-A^A„), (16) 



^AvA^;A^A7v 

2 (^AvA'„;A-,A„ - (-l)"^'' + '^"'i^-AvA'„;-A-,A„)- (17) 

From definitions (|16p . (fT7|) and relation (fT4|) the amp- 
litudes T and IJ obey the symmetry relations [23j : 

AvA^;A-yAjY I, ^} ^ — 

^(_1)AW-A«r^^_^,^^^^_^^^ (18) 



•-^ Av ; A-^ Ajv v / Av-A^; — A-yA^v 

= -(-l)'""^"t^Av-A' ;A,-A« . (19) 



For convenience, we introduce the abbreviation 

E = ^y^\' \ for the summation over the final nu- 
2 ^Ajy An 

cleon helicity indices and averaging over the initial spin 
states of the nucleon. In the following the nucleon he- 
licity indices of the amplitudes are implicit, but will 
be included when required for clarity. If Tx^,A-, appears 

without the symbol all nucleon helicity indices are 
equal to 1/2. 

For NPE amplitudes, transitions diagonal in nu- 
cleon spin (A^ = \n) are dominant. Furthermore, 
since for scattering off an unpolarized target there 
is no interference between nucleon spin-flip and non- 
spin-flip amplitudes, the fractional contribution of nu- 
cleon spin-flip NPE amplitudes to SDMEs is of the 



order of — <'/(4Af^), which is small at low t' . In this 
case, neglecting the small nucleon spin-flip amplitudes 
2^Ai/±i/2;A^=Fi/2 and using ((T8|) reduces the summation 
and averaging ^ to one term: 

Erji rji^ 1 ^ ^ rji rji^ 

JAvA^Ja'j,,A; = 2 ^ ^Av-A'„;>7>iv-'A;,-A'„:A;A„ 
Aiv A^ 

= ^Av1/2;A.,1/2Ta;,1/2;A;1/2 

+7av-1/2;A^1/27'4^-1/2:A' 1/2 



Tav 1/2; A^ 1/2^A^ 1/2; A' 1/2 



Ta^.a^T;, (20) 



where T* represents the complex conjugate quantity. 

For UPE amplitudes in general, the dominance of 
diagonal transitions (Aat = A^) cannot be proven, so 
that no relation similar to (|20|) can be derived and 

therefore ^ is always used. 

For unpolarized targets, there is no interference be- 
tween NPE and UPE amplitudes [23] as 



0, 



(21) 



following from relations ([TS]) and (jl9p without addi- 
tional assumptions. 



3.3 Spin Density Matrices of Photon and Vector 
Meson 

The photon spin density matrix normalized to unit 
flux of transverse photons cornprises the unpolarized 
([/) and polarized (L) matricecl, with Pbeam being the 
longitudinal polarization of the beam: 

^?a;a; = ^^a.a; + nea™ a;, (22) 



1 



(23) 



vMl+^e 



— i'P 



-2i<S> 



2e 



~ee'** -Ve(l + e)e 



i<l> 



(24) 







The spin density matrix pxvX{, of the produced 
vector meson is related to that of the virtual photon, 
ff^^y 1 through the von Neumann formula: 



PAvA' 



Tf ■^AvA^;A-,A„ £'a^\^ ^\*^a;,;A;A„i (25) 



2Af 



A^ A^ AjY Ajy 



^ The adjectives "(un)polarized" are used here with the 
same meaning as when applied to SDMEs. 
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where Fa^-a'„;A^a„ 



denotes the hehcity ampHtude of 
the j*N — !■ p'^N transition defined in The nor- 

maUzation factor is given by 



N = Nt + eAfh 



(26) 



with 



+ |C/iiP + |C/oi|' + |C/-iin, (27) 



2|TioP + 2|[/ion 



(28) 



Equation (|28p is obtained by using symmetry relations 
(UHl) and HI]). 

If the spin density matrix of the photon is decom- 
posed into the standard set of nine hermitian matrices 
Z"" (a = 0, 1, 8); a set of nine matrices P^i/A' 
obtained for the vector meson [23]: 



A^ A^Ajy Ajv 



y^a 771* 

A^Ajv^A^A^-f'A^A^iA^AN 



(29) 



The four matrices for a = 0, 1, 2, 3 in ((29)) describe 
vector meson production by transverse virtual pho- 
tons: unpolarized, linearly polarized in two orthogonal 
directions, and circularly polarized, respectively. For 
these cases A/'q — Nt- Vector meson production by 
longitudinal virtual photons corresponds to a = 4 in 
((29)) and Na = J^l- The interference between the amp- 
litudes of vector meson production by transverse and 
longitudinal virtual photons is described by ((29)) for 
a = 5, 6, 7, and 8 with Ma = ^/MtMl- 



3.4 Cross Sections 

The differential cross section of the reaction 7* TV 
p°iV 7T~^Tr~N is given by 



dafMW,Q^) _ fiW,Q^) 



dt d<!> d4> d cos 

X V F: 



47r 



AvA'„;A-,A„ Pa^A^I^'^) ^A 

X yiA^((/),cose) YiV.((/),cos©), (30) 



in terms of g^^y , the virtual-photon spin density ma- 
trix, the helicity amplitudes ^Av-a^:A^An describing 
the transition of the virtual photon with helicity 
to the vector meson with helicity Ay, and the spher- 
ical harmonics yim(0, cos Q), to = ±1,0 (defined as 
in pSlfHll^ ) that describe the angular distribution of 



the pions from the decay —f 7r+ +7r^. It is assumed 
here that the branching ratio of the p'^-meson decay 
into 7r+7r~ is 100%. The kinematic factor 



fiW,Q^) 



1 



(31) 



in ((30)) accounts for the fact that the flux of transverse 
photons in electroproduction is not unity (see Ref. [53] 
for the relation of the differential virtual-photon cross 
section to the differential electroproduction cross sec- 
tion). 

The singly differential cross section "^"j^"" for 
meson production is obtained by integrating ()30p over 
<P, (p, cos 6*. The integration over ^ eliminates the 
interference between contributions of transverse and 
longitudinal photons and makes the photon density 
matrix diagonal. For this case, the full differential cross 
section becomes the linear combination of the cross 
sections and °f vector meson production with 
transverse and longitudinal photons, respectively: 



da full 
dt 



do_L 
dt 



dar 



(32) 



where 



doi 
dt 



{W,Q^t) = J{W,Q^)M,{W,Q\t), (33) 



for i = L,T, where A/r and A/l are defined in ((27)) and 
((28)) . respectively. 

The "differential" longitudinal-to-transverse cross 
section ratio is defined as: 



R(W,Q^t) = -^/ 



da J ,daT Ml 



dt ' dt Mt 



(34) 



The complete representation for R in terms of helicity 
amplitudes is obtained by inserting ((28)) and ((27)) into 
((34)) . Approximate expressions for R related to SCHC 
or NPE will be discussed in section 11. 



3.5 Accessible Spin Density Matrix Elements 

For an unpolarized target and a longitudinally polar- 
ized beam, the 3-dimensional angular distribution of 
p'^ production and decay is described by 26 matrix el- 
ements Pav- A' |23) • If the experiment can be performed 
only at one beam energy, the matrix elements PavA' 



and 



cannot be disentangled, so that only 23 el- 



ements are accessible. It is customary to extract from 
the experimental data the following elements: 



p04 ^ 



(Pa.A' +ei?ptA' + 



Av AC, 



(1 + e-R) ' 

VRp' 



a = 1, 2, 3, 



(35) 



-JT+mr^ a = 5, 6, 7, 
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From now on, we will designate r^^^y and r'^^y 

(a = 1-3, 5-8) as the Spin Density Matrix Elements 
(SDMEs). 

In Appendix A, the SDMEs are ex- 

pressed in terms of NPE and UPE amplitudes, as ob- 
tained by combining (|^ and ([55)1 . 



3.6 Extraction of SDMEs from Measured Angular 
Distributions 

Measurement of the 3-dimensional p° production and 
decay angular distribution 

>V^+^(M^,Q^^,<I', 0,cos6)) 

_ da full /dof ull /or*\ 

dt d<P d(t> dcosO dt 

reveals the helicity structure of the 7* — > p'^N tran- 
sition. Its integral over the variables (j), and cos is 
equal to unity. The W, and t dependences of W^'^^ 
are contained in the corresponding dependences of 
the SDMEs r'^^^y . The full angular dependence of 

yy'^+^(<?>, 0, cos 6)), as a hnear function of the SDMEs 
'"av-a' ' ^® given in (|37ll39p as derived in Ref. [23]. (Note 
that these formulae are on the next page.) 

3.7 s-Channel Helicity Conservation. 

The measurement of SDMEs allows the determination 
of the extent to which s-channel helicity is conserved 
for a given process and kinematic conditions. SCHC 
implies that the contributions from all non-diagonal 
transitions ^Av-A^;A.^Aiv with 7^ Ay are zero. In 
terms of NPE and UPE amplitudes, only Tqo, Th, and 
Ull remain. As a consequence, all spin density matrix 
elements vanish except the unpolarized SDMEs Tqq, 
r\_i, Im{rf_j}, Re{rfQ}, Imjrj'g}, and the polarized 
ones ImjrjQ} and Re{rfQ}, as can be seen from ((771 - 
IMl) of Appendix A. If SCHC holds, SDMEs are not 
independent, as the following relations apply: 

r\_, = -Im{rti}, (40) 
Rc{r5„} = -Im{r6j, (41) 
Re{r«„} = Im{r[o}. (42) 

The measurement of SDMEs also allows the deter- 
mination of the extent to which the unnatural-parity- 
exchange mechanism is relevant for a given process and 
kinematic conditions. If natural-parity exchange dom- 
inates, so that the amplitude Un can be neglected, an 
additional relation is obtained: 

1 - rlt = '2r\_, = -2Im{rti}. (43) 



4 The HERMES Experiment 



The HERMES experiment at DESY used a 27.6 GeV 
longitudinally polarized positron or electron beam im- 
pinging on pure hydrogen or deuterium gas targets in- 
ternal to the HERA storage ring. Parts of the data set 
were collected with longitudinally or transversely po- 
larized targets, the polarization of which was flipped 
approximately every minute. The average over the tar- 
get polarization values was confirmed to be consistent 
with zero, as required for the extraction of SDMEs in 
this analysis. The lepton beam was transversely self- 
polarized by the emission of synchrotron radiation |39] . 
Longitudinal polarization at the interaction point was 
achieved by spin rotators located upstream and down- 
stream of the HERMES apparatus. For both positive 
and negative beam helicities, the beam polarization 
was continuously measured by two Compton polarime- 
ters f40[BT| . The average beam polarization for the 
hydrogen (deuterium) data set was 0.45 (0.47) after 
requiring 0.15 < Pbeam < 0.8, and the fractional sys- 
tematic uncertainty of the beam polarization was 3.4% 
(2.0%) [401I4T]. 

The data sample recorded with a longitudinally po- 
larized hydrogen (deuterium) target, representing 14% 
(38%), of the total statistics, has a residual polariza- 
tion of 0.0221 ± 0.0001 (-0.0036 ± 0.0009). The data 
sample recorded with a transversely polarized hydro- 
gen target, representing 35%, has a residual polariza- 
tion of 0.0028 ± 0.0001. The systematic uncertainty of 
the target polarization measurement is typically 0.04. 

The HERMES spectrometer is described in de- 
tail in Ref. [42]. It was a forward spectrometer in 
which both scattered lepton and produced hadrons 
were detected within an angular acceptance ±170 
mrad horizontally, and ±(40 - 140) mrad vertically. 
The scattered-lepton trigger was formed from a coin- 
cidence between three scintillator hodoscopes and a 
lead-glass calorimeter. The trigger required an energy 
of more than 3.5 GeV deposited in the calorimeter. 
The tracking system had a momentum resolution of 
~ 1.5% and an angular resolution of w 1 mrad. Lep- 
ton identification was accomplished using a lead-glass 
calorimeter, a preshower detector consisting of a scin- 
tillator hodoscope preceded by a lead sheet, and a 
transition-radiation detector. Until 1998 the particle 
identification system included a gas threshold Ccren- 
kov counter, which was replaced in 1999 with a dual- 
radiator ring- imaging Cerenkov detector (RICH) |43] . 
Combining the responses of these detectors in a likeli- 
hood method leads to an average lepton identification 
efhciency of 98% with a hadron contamination of less 
than 1%. 
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W"+^{'P,(t>,cos0) = {<P, (j), cos 0) {$,(!), cos 0), 



(37) 



W^('5,0,cos©) = hl-r'ot) + hsr°t>-l)cos^0~V2Re{r1t}sm20cos- 



— e cos 2<?>^r}i sin^ + t-qq cos^ — \/2Re{riQ} sin 26* cos cji 

— e sin 2<?>^\/2Ini{rio} sin 20 sin + lia{ri_i} sin^ 6* sin 2<f) 



r?iiSin^ecos 2(^ 
2<^) 



1 • 2 ^ 

- sm 6* cos 



+ ^/Ml+7jcos<p( rfi sin^ + roo cos^ 6* — V2Re{rfo} sin 29 cos (/> — ri_i sin 0cos2(j) 



yV''{$,4>,COS0) = -^Pbeam 



+ V2e(l + e) sin *(^\/2Im{r?o} sin 20 sin (/i + Im{r-i_i} sin^ sin 2(/) 

\/l — ^V^Ini{rio} sin 20 sin + Im{r-i_i} sin^ sin 2(/>^ 
+ V2e(l-e) cos<?(^^I: 
+ A/2i(r^sin<?(^r?isin^6> + r'* 



[{rJo} sin 20 sin + Im{ri_i} sin^ sin 2(t> 



r-Qo cos^ — \/2Re{rio} sin 20 cos — rf_i sin'^ 



cos 20) 



(38) 



(39) 



5 Data Analysis 

5.1 Exclusive p° Events 

Events accepted for the analysis are required to fulfill 
the following criteria (see Refs. [T2[|44] for details): 

— three tracks originate from the target and are 
recorded in the spectrometer; 

— two oppositely charged hadrons and one lepton 
with the same charge as the beam are identified 
through the likelihood analysis of the combined 
responses of the four particle-identification detec- 
tors [42] : 

— the reconstructed virtual photon satisfies the kine- 
matic constraint 1 GeV^ < < 7 GeV^; 

— the p° meson is selected by requirements on the in- 
variant mass of the two hadrons of opposite charge: 
0.6 GeV < M^+^- < 1.0 GeV when both hadrons 
are assumed to be pions, and the veto constraint 
Mx+K- > 1-06 GeV, the latter under the hypoth- 
esis that both hadrons are kaons. The veto con- 
straint excludes contamination from —>■ 
decay. Two-pion invariant mass distributions in the 
HERMES acceptance for proton and deuteron data 
are presented in Fig. IH A detailed description of 
the invariant mass peak of exclusive p° events is 
published in Refs. [T51IT8] and also presented in 
Ref.|12j. The distribution of these events in both 
AE and t' is presented in Fig. [5l 

— exclusive events are selected by the requirement 
-1.0 GeV < AE < 0.6 GeV (called the "exclusive 
region" in the remainder of the text). The applica- 
bility of such a constraint was explained in detail 
in Ref. [HI, as weU as in Refs. [IMHEU • In the 



AE spectrum the resolution due to instrumental 
effects ranges between 0.26 and 0.38 GeV depend- 
ing on the spectrometer configuration. 
— the "final event sample" of p° events is obtained 
from the sample of exclusive events by the addi- 
tional requirement —t' < 0.4 GeV'^. This require- 
ment ensures that the semi-inclusive background 
does not exceed a level of about 10% in the kine- 
matic bins of and t' presented below. 

After the application of all the above requirements, 
the entire kinematic region contains 16362 p*^ events 
from hydrogen and 25940 events from deuterium, which 
are used in the subsequent physics analysis. 

5.2 Backgrounds for Exclusive p° Events 

In exclusive vector meson production, the target nu- 
cleon remains intact. At HERMES the recoiling tar- 
get nucleon was not detected and hence, given the ex- 
perimental resolution, also a certain number of non- 
exclusive events will satisfy the requirements for ex- 
clusive events. They appear as background remaining 
underneath the AE peak. 

a) Background from Semi-Inclusive Deep-Inelastic Scat- 
tering 

Background events originate mainly from fragmenta- 
tion processes in Semi-Inclusive Deep-Inelastic Scat- 
tering (SIDIS), in which the final state contains a pair 
of oppositely charged hadrons in the spectrometer. 
Only a small fraction of this background passes the 
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Fig. 4. Two-pion invariant mass distribution in the spec- 
trometer acceptance, fitted with a skewed Breit-Wigner 
function (solid line). The dotted line represents a constant 
background contribution. 




Fig. 5. The distribution of mesons in AE, and t' from 
the 1996-2005 hydrogen data sample in the range 1 GeV^ < 
< 7 GeV^ 



above-described AE and —t' requirements for exclu- 
sive p° production. 

The amount of SIDIS background and its angular 
distributions in the exclusive region can not be de- 
termined with the present apparatus. Therefore, the 
PYTHIA code gS] tuned for HERMES kinematic con- 
ditions [46 l l47 l l48] is used. Exclusive processes were 
excluded from the simulated sample. The simulated 
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Fig. 6. Distribution of AE for the 1996-2005 hydrogen 
data sample shown for four intervals in i! after application 
of all event selection requirements, except the one for AE. 
The shaded areas represent the SIDIS background calcu- 
lated by PYTHIA, normalized to the data in the region 
2 GeV < AE < 20 GeV. 



SIDIS events were passed through the same chain of 
kinematic requirements as the experimental ones. Very 
good agreement between the experimental data and 
the simulation is observed for the shape of the AE 
distributions in the region AE > 2 GeV for each kine- 
matic interval in (5^, x or t', and for both targets. 
This agreement in shape is demonstrated in Fig. [6l 
which shows four intervals in t' as an example. Since 
no absolute normalization between data and simula- 
tion is required to determine the SDMEs (as shown 
in the next section), for every kinematic interval the 
fractional background contribution ftg in the exclusive 
region can be obtained by normalizing the simulation 
to the data in the region 2 GeV < AE < 20 GeV. 
That is, 



(44) 



where N^^%, Ni% 



and N' 



data 



e.ci > N^^ci are the total 
number of measured and simulated events at 2 GeV < 
AE < 20 GeV and in the exclusive region, respec- 
tively. This contribution amounts to 8% for the entire 
kinematic region and ranges between 3 and 12% in the 
different kinematic intervals. These values will be used 
for the subtraction of the SIDIS background angular 
distributions, as described in section 6.2. 



b) Background from Non-resonant Exclusive tt+tt Pairs 

The contribution of non-resonant tt+tt^ production 
and its interference |49| with resonant p'^ — > tt"'" tt" 
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production are determined using a modified Breit- 
Wigncr fit to the invariant mass distribution. We do 
not distinguish between resonant and non-resonant 
contributions in production, following the prac- 
tice of previous experimental [9j[50] and theoretical 
publications P^[5Tl[551F5i] . Therefore the present data 
were not corrected for non-resonant background, which 
amounts to ^ 4 - 8% depending on the kinematics [151 
[T7HT2] . Note that the contribution of p'^-uj interference 
has been found to be negligible [T5] . 



c) Background from Proton-Dissociative Processes 

Another possible background consists of events in 
which the target proton is excited to some other bary- 
onic resonance, which then decays to a proton and 
typically a soft pion. In the absence of a recoil de- 
tector, such events cannot be distinguished, but their 
contribution to the exclusive p'^ production cross sec- 
tion at HERMES was found to be small (4± 2%) [15] . 
No correction has been applied for this background, as 
the extracted SDMEs were found to change negligibly 
at a value of AE=Q.2 GeV, where this background is 
strongly suppressed. This approach is supported by re- 
sults from ZEUS, where the decay angle distributions 
of the proton-dissociative reaction have been found to 
be consistent with those of exclusive events [9lfT0]. 



6 Extraction of SDMEs 

6.1 Maximum Likelihood Method 

In each bin of or t' , or the entire acceptance, 
the SDMEs are obtained by minimizing the differ- 
ence between the 3-dimensional (cos 0, (/), produc- 
tion and decay angle distribution of the experimen- 
tal events and that of a sample of fully reconstructed 
Monte Carlo events, using a binned maximum likeli- 
hood method. For the Monte Carlo simulation, events 
are generated isotropically in (^, 0, cos 0) using the 
rhoMC generator [T5 l fT6 l l44] for exclusive p° produc- 
tion, simulated in the instrumental context of the spec- 
trometer, and passed through the same reconstruc- 
tion chain as the experimental data. Introduction of 
estimated misalignments of the spectrometer into the 
Monte Carlo simulation used for the SDME extraction 
was found to have a negligible effect on the results. The 
variables cos 0, (f), and arc binned in 8 x 8 x 8 cells. 
The content of each of the 512 cells is weighted using 
([37]) . whereby the 23 matrix elements are treated as 
free parameters. The number of events di in each cell 
is assumed to obey a Poisson distribution: 



with mean value cniti^, where cn = "^'j) 
is a normalization factor accounting for the difference 
in the total number of events in the data {dj ) and simu- 
lated (TOj ) sample, and m[ is the (re) weighted number 
of simulated events in cell i. The likelihood function is 
then defined as [5T] 



cells 



LW - Y[Pid^,CN{^)mr{X)) 



(46) 



where A represents the 23 fit parameters that are the 
23 SMDEs. The best fit parameters were determined 
by maximizing the logarithm of the likelihood func- 
tion, 

InL(A) = [^'^ 1^ (cw(A)m^(A)) - CAr(A)m^] 

i 

+ constant, (47) 

or equivalently by minimizing — ln_L(A). 

The minimization itself and the uncertainty calcu- 
lation are performed using the MINUIT package [52] . 
In the fitting procedure the samples with negative and 
positive beam helicity are fitted simultaneously. The 
values of Pcr degree of freedom {x^ /d.o.f.) for the 
16 kinematic intervals {Q^, t' or x), calculated af- 
ter completing their likelihood fits, range between 0.6 
and 1.2 for 8 x 8 x 8 — 23 degrees of freedom. For ev- 
ery SDME, the averages of the SDMEs extracted from 
the two separate beam helicity samples are found to 
be consistent with each other and the result from the 
common fit. 

In Fig. [7] one-dimensional angular distributions are 
shown for the hydrogen and deuterium data samples, 
where the positive-helicity sample is chosen as repre- 
sentative. In addition to distributions in cos 0, (j), and 
^, the angular distribution in ■(/; = — is shown, 
which embodies the entire azimuthal dependence in 
the case of SCHC. For each panel, the data are com- 
pared with the isotropic input distributions as modi- 
fied by instrumental effects such as acceptance, track- 
ing resolution, and reconstruction efficiencies, as well 
as the one-dimensional projections of the fitted 3-di- 
mensional angular distribution. These projections are 
clearly in agreement with the data. 

6.2 Background Subtraction 

Before fitting SDMEs to the (cos 0, 0, ^) angular dis- 
tributions, the SIDIS background in the exclusive re- 
gion is subtracted. This subtraction is performed sep- 
arately for each interval in t' , or x. In a given 
(cos©, </),<?) cell, the number of background events in 
the exclusive region is calculated as follows: 



P{di,CNm[) 



{cN■m'^ 



id. 



d^l 



(45) 



j\Tdata 

IVfbg ^ ivrMC ^"2-20 
cell ^^cell j^MC ' 



(48) 
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Fig. 7. Angular distributions with common arbitrary normalization for p meson production and decay. Data points 
represent the positive helicity sample of the proton (deuteron) data in the left (right) half of the figure. The dotted lines 
represent isotropic input Monte Carlo distributions as modified by the HERMES acceptance, while the dashed lines are 
the results of the 23-parameter fit. Here 16 bins are chosen for a more detailed comparison. The data correspond to the 
full kinematic region of the analysis. 



where the number of SIDIS Monte Carlo events in a 
given cell is iV*/,f , while N^lfo and Nl% are defined 
as in equation 



6.3 Systematic Uncertainties 

a) Background Subtraction 

The systematic uncertainty of the background subtrac- 
tion procedure is estimated as the difference between 
the SDMEs obtained with and without any background 
correction. 



b) rhoMC Input Parameters 

The SDME extraction procedure starts from isotropic 
distributions in cos0, (p, and generated by rhoMC, 
as explained above. The parameterization of the to- 
tal electroproduction cross section in rhoMC is chosen 
in the context of a VMD model that incorporates a 
propagator-type dependence, and also contains a 
dependence on R{Q^). As the HERMES spectrometer 
acceptance depends on , different input parameters 
result in sUghtly different reconstructed angular dis- 
tributions. The corresponding systematic uncertainty 
of the resulting SDMEs is obtained by varying these 
parameters within one standard deviation in the total 
uncertainty of the parameters given in Refs. [151ll6j . 



The total systematic uncertainty is obtained by 
adding in quadrature the uncertainty from the back- 
ground subtraction procedure and that due to the un- 
certainty in the rhoMC input parameters, which are 
approximately of equal size. 



c) Further Systematic Studies 

Several further studies using generated and reconstruc- 
ted event samples arc performed to estimate possible 
systematic uncertainties: 

i) A consistency check of the extraction method is 
performed by using several known sets of SDMEs as 
input to the rhoMC [TMTMH] simulation and com- 
paring the SDMEs extracted from the simulated data 
with those used as input to the rhoMC generator. 
First, isotropic angular distributions were simulated, 
corresponding to all SDMEs vanishing except Tqq = i . 
Alternatively, events were generated assuming SCHC, 
implying that only five unpolarized and two polarized 
SDMEs are non-zero, as explained in section 3.7. Fi- 
nally, the extracted 23 SDMEs arc used as input pa- 
rameters. In all cases, good agreement is found be- 
tween input and extracted SDMEs at the given level 
of statistical accuracy. 

ii) Several tests are performed to ensure that the 
choice of the {cos0, (fijip) cell size does not bias the 
results of the maximum likelihood procedure. A sam- 
ple of about 40000 simulated events with angular de- 
pendences determined by the (normally) extracted 23 
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SDMEs is fitted after binning the data in several num- 
bers of angular cells, varying from 5 x 5 x 5 to 12 x 12 x 
12. The x^/d.o.f. calculated between sets of SDMEs 
extracted with 8x8x8 and 12x12x12 binning is 
0.14. Hence the cell size used in the maximum likeli- 
hood procedure is not treated as a source of systematic 
uncertainty. 

iii) Variations of the restrictions on Af^+jr- , t', and 
AE result in slightly different amounts of SIDIS back- 
ground. The resulting systematic uncertainty is much 
smaller than that estimated for the background sub- 
traction procedure, and hence is neglected. 

iv) In the SDME extraction procedure, only the 
shape of the 3-dimensional angular distribution mat- 
ters. As events in which a radiative photon is emit- 
ted with an energy larger than 0.6 GeV are removed 
from the analysis by the constraint AE < 0.6 GeV, 
the impact of radiative effects on the shape of the 3- 
dimensional angular distribution is strongly reduced. 
Two approaches are used to quantify this effect. First, 
the DIFFRAD code is used to calculate the radiative 
effects for exclusive p° production [53ll54] , as was done 
also in Refs. pilT^ . As the emission of a real photon by 
the positron alters the direction of the virtual photon, 
the angle between lepton scattering plane and p° 
production plane also changes. The effect of a small 
variation (< 2.5 %, as suggested in Ref. [54]) of the 
shape of the ^-distribution is studied by re-weighting 
the isotropic input angular distribution. The differ- 
ence between SDMEs obtained with and without re- 
weighting is found to be less than 0.0012 for all SDMEs 
(x^ /d.o.f. < 0.1), i.e., radiative effects are negligible. 

As an independent cross check, radiative effects on 
the extracted SDMEs are studied using a Monte Carlo 
simulation of exclusive production with events from 
the PYTHIA generator [45]. Two statistically inde- 
pendent isotropic angular distributions are generated, 
with and without the emission of radiative photons. A 
set of SDMEs is extracted from the (real) data sample 
for each of these isotropic input angular distributions. 
The difference between the resulting SDMEs is statis- 
tically indistinguishable (x^/d.o.f. < 0.2). 

As a further check we use the extracted 23 SDMEs 
as input parameters to rhoMC and compare the shapes 
of the simulated distributions with the data. In order 
to restrict the comparison to exclusive p*^ events, prop- 
erly normalized SIDIS background distributions from 
PYTHIA are subtracted from the data. In the maxi- 
mum likelihood fit method, the extraction of SDMEs 
only requires simulated event distributions normalized 
to the data. The shape comparison reveals good agree- 
ment for all variables, some of which are presented in 
Fig.[Hl 
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Fig. 8. Shape comparison of the distributions in W , Q , 
t' and , the momentum of the vr^ from decay in the 
laboratory system, for the hydrogen data sample (squares) . 
The shaded areas show rhoMC results using the extracted 
23 SDMEs as input for the simulation, normalized to the 
data. Background has been subtracted from the data. 

7 Results on SDMEs Integrated over the 
Entire Kinematic Region 

7.1 Classification of SDMEs 

The presentation of the extracted SDMEs will be based 
on a hierarchy of NPE helicity amplitudes: 



iTn, 



iTiil » iToil > iTiol > |ri_ 



(49) 



A similar hierarchy was discussed for the first time in 
Ref. [2^ . In perturbative QCD this is valid at asymp- 
totically high Q^. It was experimentally confirmed at 
the HERA collider [9ti 10^11] and discussed in Refs. [M] 
[3T1I33] . In the following it will be shown that it applies 
also at values typical of the HERMES experiment. 

The extracted 23 SDMEs are categorized into five 
classes according to the hierarchy shown above. Class 
A comprises SDMEs dominated by the helicity-conser- 
ving amplitudes Tqq and Tn which describe the tran- 
sitions 72 p'j^ and -f^ p^, respectively. Class B 
contains SDMEs that correspond to the interference of 
the above two amplitudes. Class C consists of all those 
SDMEs in which the main term contains a contribu- 
tion linear in the s-channel helicity non-conserving am- 
plitude Toil corresponding to the 7^, — > transition, 
except for a term involving Tqq for which the Tqi contri- 
bution is quadratic. The classes D and E are composed 
of the SDMEs in which the main terms contain a con- 
tribution linear in the small helicity-flip amplitudes 
^10 ill ^ Pt) ^^^d Ti_i (7lj, -> p!^), respectively. 
Equations (|77ll99p in Appendix A show the represen- 
tation of all the SDMEs in terms of helicity amplitudes 
ordered according to the five classes defined above. 
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Table 1. Values of p SDMEs for proton and deuteron data ordered in classes by horizontal lines according to the 
expected hierarchy of helicity amplitudes. Elements with a = 3, 7, 8 are polarized. The first uncertainty is statistical, 
the second systematic. The statistical significance of the magnitude of each SDME expected to vanish in the case of 
SCHC is shown in column three or five as its absolute value in units of standard deviations of its total uncertainty, 
denoted as SDME/fot. Similarly, the four bottom rows of the table show deviations from zero of certain combinations 
of SDMEs that would thereby violate SCHC and NPE dominance (see text). 





proton 


deuteron 


element 


SDME ±stat ± syst 


SDME/tof 


SDME ±stat ± syst 


SDME/tof 




0.412 ± 0.010 ± 0.010 




0.416 ± 0.007 ± 0.013 






0.246 ± 0.011 ± 0.019 




0.247 ± 0.008 ± 0.014 




Im ri_i 


-0.227 ± 0.010 ± 0.012 




-0.234 ± 0.008 ± 0.019 




Re rfo 


0.161 ± 0.004 ± 0.010 




0.165 ± 0.003 ± 0.010 




Im rfo 


-0.167 ± 0.003 ± 0.009 




-0.156 ± 0.003 ± 0.006 




Im r^o 


0.112 ± 0.021 ± 0.011 




0.104 ± 0.016 ± 0.004 




Re rfo 


0.074 ± 0.019 ± 0.006 




0.114 ± 0.014 ± 0.009 




Re rlt 


0.031 ± 0.004 ± 0.008 


3.5 


0.030 ± 0.003 ± 0.008 


3.5 


Re rJo 


-0.032 ± 0.007 ± 0.012 


2.3 


-0.020 ± 0.006 ± 0.010 


1.7 


Im rfo 


0.022 ± 0.007 ± 0.015 


1.3 


0.014 ± 0.006 ± 0.017 


0.8 




0.109 ± 0.009 ± 0.009 


8.7 


0.111 ± 0.007 ± 0.008 


10.4 


roo 


0.011 ± 0.019 ± 0.008 


0.6 


-0.038 ± 0.015 ± 0.016 


1.7 


Im rfo 


-0.017 ± 0.015 ± 0.004 


1.1 


0.031 ± 0.011 ± 0.003 


2.7 


roo 


0.035 ± 0.050 ± 0.010 


0.7 


0.053 ± 0.038 ± 0.006 


1.4 


rfi 


-0.01 ± 0.003 ± 0.013 


1.2 


-0.021 ± 0.003 ± 0.013 


1.6 


rf-i 


0.005 ± 0.004 ± 0.006 


0.7 


0.013 ± 0.003 ± 0.006 


1.9 


Im rf_i 


-0.002 ± 0.004 ± 0.007 


0.3 


-0.007 ± 0.003 ± 0.006 


1.0 


Im rf_i 


-0.035 ± 0.030 ± 0.005 


1.2 


-0.058 ± 0.023 ± 0.009 


2.3 


rii 


0.036 ± 0.024 ± 0.001 


1.6 


0.026 ± 0.018 ± 0.003 


1.4 


ri-i 


0.019 ± 0.029 ± 0.005 


0.6 


-0.066 ± 0.022 ± 0.008 


2.8 


ri-i 


-0.011 ± 0.005 ± 0.005 


1.5 


-0.002 ± 0.004 ± 0.008 


0.2 


rii 


-0.025 ± 0.007 ± 0.008 


2.3 


-0.002 ± 0.006 ± 0.013 


0.1 


Im rf_i 


-0.024 ± 0.018 ± 0.001 


1.3 


-0.004 ± 0.014 ± 0.004 


0.3 


relation 


SCHC? 


SCHC? 


r}_i + Im rf_i 


0.018 ± 0.012 ± 0.011 


1.1 


0.013 ± 0.009 ± 0.014 


0.8 


Re rfo+ Im rfo 


-0.006 ± 0.004 ± 0.001 


1.5 


0.010 ± 0.003 ± 0.005 


1.7 


Im rfo- Re rfo 


0.038 ± 0.029 ± 0.006 


1.3 


-0.011 ± 0.022 ± 0.006 


0.5 


relation 


SCHC and NPE? 


SCHC and NPE? 


1 9r-l 

J- — '^OO ^ ^ri_i 


0.097 ± 0.017 ± 0.046 


2.0 


0.091 ± 0.013 ± 0.038 


2.3 



7.2 Representation of the Integrated Data 

Separate maximum likelihood fits to the proton and 
deuteron data samples are performed in the entire kine- 
matic region: 1 GeV^ < < 1 GeV^, 3 GeV < 
W < 6.3 GeV (corresponding to 0.03 < xb < 0.25), 
and GeV^ < ~t' < 0.4 GcV^. The resulting p° me- 
son SDMEs r'^^y are listed in Tab. [T]and displayed 

in Fig. [9l ordered according to the classes described 
above. The statistical uncertainties are larger for the 
eight polarized SDMEs (shown in the shaded areas of 
the figure) due to the non-unity of the beam polar- 
ization and the kinematic suppression factor \/l — e 



(see (|39| ). In order to facilitate comparison with a re- 
cently introduced new representation of SDMEs pS] , 
the proton SDMEs in that representation are shown 
in Tab. [H of Appendix E. 

In Fig. [9] the SDMEs are shown multiplied by cer- 
tain factors in order to allow their comparison at the 
level of dominant amplitudes (see ()77II99|) ). For all 
classes numerical factors are chosen in such a way 
that the coefiicient of the dominant terms is equal 
to unity. The plotted representatives for the elements 
of class A are chosen so that their main terms are 
equal to jTiip/A/"; in particular this requires that the 
term 1 — rSn be chosen. The coefficients for class B 
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Fig. 9. The 23 SDMEs extracted from p" data: proton (squares) and deuteron (circles) in the entire HERMES kinematics 
with {x) = 0.08, (Q^) = 1.95 GeV^ {-t') = 0.13 GeV^. The SDMEs are muhiplied by prefactors in order to represent 
the normahzed leading contribution of the corresponding amplitude (see (|77l99p ). The inner error bars represent the 
statistical uncertainties, while the outer ones indicate the statistical and systematic uncertainties added in quadrature. 
SDMEs measured with unpolarized (polarized) beam are displayed in the unshaded (shaded) areas. The vertical dashed 
line at zero is indicated for SDMEs expected to be zero under the hypothesis of SCHC. 



are chosen to have the main contribution to the plot- 
ted representatives for the unpolarized and polarized 
SDMEs equal to Rc{TnT^o}/-^ and Im{TiiT*o}/7V, 
respectively. This corresponds to the general rule that 
is applicable to classes B to E: the dominant contribu- 
tion of the unpolarized (polarized) element presented 
in Fig. [5] is proportional to the real (imaginary) part 
of a product of two amplitudes. Class C contains the 
main terms TqiT^q/M (for r^o/^/2 and r§o/\/2) and 



TaiTii/Af. The dominant contributions for classes D 
and E contain terms TioT^^/Af and Ti-iT*^/J\f, re- 
spectively. 

Given the scaled SDMEs in Fig. [9l it easily can 
be seen that the two unpolarized SDMEs of class B 
have large values, similar to those of class A. This 
suggests the presence of a substantial interference be- 
tween the two dominant amplitudes Tqq and Tn. The 
two polarized class B SDMEs are significantly non- 
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zero for proton and deuteron as well. It is also seen 
that the values of elements in class C that contain 
the dominant term TqiT*i are similar for the unpo- 
larized SDMEs (Re{r?;J}, Re{rio}, Im{r?o}). Those 
unpolarized class C elements measured with good ac- 
curacy, Rejrj'Q} and Vqq, are much smaller than the 
class B SDMEs, whereas the unpolarized class C ele- 
ments are larger than the unpolarized class D and class 
E SDMEs. This shows that the anticipated hierarchy 
is supported by the data. For class D SDMEs, slightly 
positive (negative) values are observed in the polarized 
(unpolarized) case. Finally, values of class E SDMEs 
for the proton target tend to deviate from zero, while 
those for the deuteron ones are consistent with zero. 

We note that no significant difference is found be- 
tween the sets of SDMEs for proton and deuteron, 
as a x^/d.o.f. — 22.5/23 is obtained taking into ac- 
count the total imcertainties. Hence there appears to 
be no indication of significant contributions of sec- 
ondary reggeons with isospin 1=1 and natural parity. 



7.3 Test of the SCHC Hypothesis 

As explained in section 3.7, only the following seven 
SDMEs are not restricted to be zero in the case of 
s-channel helicity conservation: Tqq, rl_i, Im{r^_]^}, 
Re{rfo}, lm{r%} and Im{r[o}, Re{r^o}. All other 
SDMEs are required by SCHC to be zero. The mag- 
nitudes of their measured offsets from zero, expressed 
in units of the standard deviation of their uncertainty, 
are shown in one of two separate columns of Tab. [1] 
next to the respective SDME. Several elements are in- 
consistent with the hypothesis of SCHC, in particular 
several members of class C. 

The SDME Tqq is observed to be non-zero at the 
level of nine (ten) standard deviations in the total 
uncertainty for the proton (deuteron) result, proving 
s-channel helicity non-conservation. This was already 
observed earlier by the HERA collider experiments [9l 
[TT] at a lower significance level, and with high sig- 
nificance very recently |TU]. For the first time, HER- 
MES observes s-channel helicity non-conservation also 
in other class C SDMEs, in particular Re{r°Q}. 

The polarized elements Tqq and Imlrfg}, related 
to the terms Im{ToirQ*Q} and ImjToirf]^} respectively 
(|89l90p . are extracted using a longitudinally polarized 
lepton beam for the first time. Both elements are con- 
sistent with zero (Figs. [9l [131) within the uncertainties. 

The relations imposed by the hypothesis of SCHC 
(|40M2[) are satisfied within about one standard devia- 
tion of the total uncertainty, as can be seen from the 
corresponding rows of Tab. [TJ The sensitivity of these 
relations to SCHC is low. In the case of the relation 
((ID)) only the contributions of small double-helicity-flip 
amplitudes (see (|78l79p ) violate SCHC. For the rela- 
tions (|41ll42p . equations (|80ll83p show that the largest 
SCHC amplitude Tqq is multiplied by the smallest Ti_i 



amplitude in the terms that violate SCHC. The re- 
lation corresponding to the combined hypotheses of 
SCHC and NPE dominance (j43|) is marginally vio- 
lated by two standard deviations in the total uncer- 
tainty. In evaluating the uncertainties of these rela- 
tions, correlations between the corresponding elements 
(see Tabs. [T5ll6p . are taken into account. 



7.4 Phase Difference between Tn and T 



00 



The phase difference 6 between the amplitudes Tn and 
Too can be evaluated as follows: 



cos 5 



2V?(Re{rfo} - lm{r%}) 



,04 



(1- 



'■l-i - M'-Li}) 



(50) 



This results in \6\ = 26.4± 2.3stQt ±4.9syst degrees for 
the proton and \6\ = 29.3± 1.6stat±3.6syst degrees for 
the deuteron (see Fig. [T^ . Using polarized SDMEs, 
also the sign of 5 can be determined: 



3in 5 



2Vi(Re{rfo}+Im{rIo}) 



rj_i - Im{r?_i}) 



(51) 



Equations ([50)) and ([51)) are derived in Appendix B. 
Second order contributions of spin-flip amplitudes are 
neglected in obtaining these formulae. 

Using ([5T)) it is determined, for the first time, that 
the sign of 5 is positive: 5 = 
degrees for the proton and 5 



30.6 ± 5.0,t„t ± 2.4. 



36.3 ±3.9, 



±1.7 



syst 
syst 



for the deuteron. These values are consistent with each 
other and their magnitudes are in agreement with the 
results obtained with ([50)) for cos (5. 

We note that in the GPD-based model of Ref . [35] , 
the phase difference between the amplitudes Tn and 
Too is found to have a value of about 3 degrees. This 
appears to be inconsistent with the HERMES results 
and also, to a lesser extent, with the HI results [TT) : 
the two experimental results agree within their total 
uncertainties. 



8 The and t' Dependences of the 
SDMEs 

In the following, the dependences are presented in 
four bins, where the first bin is defined by 0.5 GeV^ < 
g2 < 1 GeV^ with (Q2) ^ 0.83 GeV^. For the t' depen- 
dences, also shown in four bins, only data with > 1 
GeV^ are included. The average value of t' is almost 
independent of and W. 



8.1 Class A: Dominant Transitions 7£ and 



7t 



Class A comprises SDMEs corresponding to the dom- 
inant transitions, — > pj^ and 7^ — > p^, described 
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Fig. 10. and t' dependences of class A SDMEs describing the dominant transitions 7£ — > p% and — » p^. Filled 
squares (circles) correspond to proton (deuteron) data. Total uncertainties are depicted, calculated as statistical and 
systematic uncertainties combined in quadrature. Deuterium data points are presented with a small horizontal offset 
here and in Figs. pi)ll8|) to improve their visibility. 
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Fig. 11. and t' dependences of the class B SDMEs describing the interference of the dominant transitions 7I, — > pi 
and ■Jt ^ Pt- Filled squares (circles) correspond to proton (deuteron) data. Total uncertainties are depicted, calculated 
as statistical and systematic uncertainties combined in quadrature. 
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Fig. 12. The Cf' dependence of the phase difference & be- 
tween Til and Too amphtudes calculated according to (|50p 
for the proton (filled squares) and deuteron (filled circles) 
data. The values of 5, for yields integrated over the range 
1 GeV^ < < 7 GeV^ are shown as open symbols. The 
inner (outer) bars represent the statistical (total) uncer- 
tainty. 



An indication of a Q dependence of the phase 
difference 8 between the amplitudes T\\ and Too (see 
([50]) ) is presented in Fig. [121 The result of a fit with a 



linear (5^ dependence has a x^/rf-o./. = 1.41/2 (1.42/2) 
for the proton (deuteron) data, which is smaller than 
the fit result with no dependence: ^/d.o.f. = 
4.52/3 (4.38/3). Note that at the lowest Q^, the value 
of 5 has the largest systematic uncertainty due to the 
rapidly falling acceptance of the HERMES spectrome- 
ter. No t' dependence of 5 is observed, for either target. 



8.3 Class C: Helicity-Flip Transition 7^ 



Class C consists of all those SDMEs with the main 
term containing a product of the s-channel helicity vi- 
olating amplitude Tgi describing the helicity-flip tran- 
sition 7y —* p^, and the complex conjugate of Too, 
Til or Tqi, (see (|84j|90p ). No clear dependence is 
observed for class C SDMEs (see Fig. [HI). 

As suggested by general properties of helicity-flip 
kinematics at low t' values, a dependence on \J —t' that 
monotonically increases from zero is expected for the 
amplitude Toi |25| . This is consistent with the mea- 
sured SDMEs containing this amplitude, as is clearly 
seen for rgo and Re{rio} in the third row of Fig. [131 



by the amplitudes Tn and Too. The and t' de- 
pendences for the class A SDMEs 1 — rpQ, 7'i_i, and 
Im{ri_i} are shown in Fig. [TOl The three elements 
exhibit somewhat similar dependences. They are 
found to be approximately constant over the measured 
t' range, as also observed by ZEUS [10] for at av- 
erage values of 3 and 10 GeV^. Such a t' indepen- 
dence indicates similar f'-slopes for longitudinal and 
transverse components of the vector meson production 
cross section. 

We note that there is good agreement between the 
highest points of the HERMES proton data and 
the GPD-based model calculations of Ref. [3S] . 



8.4 Class D: Helicity-Flip Transition 7^ p%. 

Class D consists of SDMEs for which the main terms 
in (|9m96p contain a product of the small helicity-flip 
amplitude Tio with the complex conjugate of T\\. Un- 
polarized (polarized) SDMEs represent the real (imag- 
inary) part of this product. Correspondingly, they de- 
scribe the interference of the helicity-flip transition 
72 — > with the helicity-conserving transition 7^ — *■ 
p\. Figure [Ml shows that the class D SDMEs depends 
only weakly on and t\ and are consistent with zero 
as —t' 0, as expected. 



8.2 Class B: Interference of 7^ ^""^ 7t ^ Pt 

Transitions 

Class B comprises SDMEs describing the interference 
of the dominant transitions 7^ — > p\ and 7^ — > 
i.e., those corresponding to a product of the amplitude 
Til and the complex conjugate of Too. Polarized (un- 
polarized) SDMEs correspond to the real (imaginary) 
part of this product. 

Figure [TT] shows the and t' dependence of these 
SDMEs. It is apparent that the SCHC relations ([41]) 
and (|42p are approximately fulfilled over the mea- 
sured kinematic ranges. Considering (|80II83|) . this im- 
plies that contributions of single- and double-helicity- 
flip amplitudes are small. 



8.5 Class E: Double Helicity-Flip Transition 

* 
7-T Pt 

Class E consists of the SDMEs with the main term 
describing the interference of the transition 7!^, p^ 
with 7y — > This corresponds to a product of the 
double spin-flip amplitude Ti_i with the complex con- 
jugate of the helicity-conserving amplitude Tn. Unpo- 
larized (polarized) SDMEs represent the real (imag- 
inary) part of this product. Their and t' depen- 
dences are presented in Fig. 1151 where it can be seen 
that the class D SDMEs depend only weakly on 
and t', and are consistent with zero as —t' 0, as 
expected. 



19 





Fig. 13. and t' dependence of the class C SDMEs describing tlie interference of tiie lielicity-fiip transition 7^ — > p'i 
and one of the dominant helicity-conserving transition. Filled squares (circles) correspond to proton (deuteron) data. 
Total uncertainties are depicted, calculated as statistical and systematic uncertainties combined in quadrature. 
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Fig. 14. and t' dependences of class D SDMEs describing the interference of the lielicity-fiip transition j1 — > 
and the dominant transition Jt ^ Pt- Filled squares (circles) correspond to proton (deuteron) data. Total uncertainties 
are depicted, calculated as statistical and systematic uncertainties combined in quadrature. 
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Fig. 15. and t' dependences of the class E SDMEs describing the interference of the double-helicity-fiip transition 
7*y and the dominant transition 7^ p^. Filled squares (circles) correspond to proton (deuteron) data. Total 

uncertainties are depicted, calculated as statistical and systematic uncertainties combined in quadrature. 



9 Unnatural-Parity Exchange 

For production on the proton, or incoherent pro- 
duction on the deuteron, the existence of unnatural- 
parity exchange can be tested by evaluating the fol- 
lowing combinations of SDMEs: 

u, = l^rlt + 2rlU - 2ri\ - 2r\_,, (52) 

U2 = rJi+rf_i, (53) 
"3 = r?i+rti. (54) 

If UPE is absent, all three combinations are expected 
to vanish without resort to SCHC. A non-zero result 
for 

ui = ^{4e|i7ioP + 2\Uii + C/„iip}/AA (55) 

indicates the existence of UPE contributions, while the 
value for 

U2 + m3 = \/2^{(;7ii + C/_ii)*;7io}/AA (56) 

can vanish despite the existence of UPE contributions. 
Such a behavior can be explained if a hierarchy exists 
also for unnatural-parity-exchange amplitudes: 



This hierarchy is analogous to ((49|) and can be assumed 
to be a general property of UPE amplitudes. 

The proton result m = 0.125±0.021stQt±0.0505yst 
for the entire HERMES kinematic region differs from 
zero at a level of more than two standard deviations 
in the total uncertainty, suggesting the existence of 
unnatural-parity-exchange contributions. The deuteron 
result ui — 0.091 ± Q.QlQgtat ± 0.046syst also exceeds 
zero, but with smaller significance. Note that for both 
targets, systematic uncertainties strongly dominate. 
The dependences on and t' of ui for the proton 
and deuteron are presented in Fig. [16] and Tab. [TOl 
Although the uncertainties are large due to the large 
number of SDMEs involved in relation ([5^ . all mea- 
sured values of ui arc positive over the whole kinematic 
range. For the calculation of uncertainties in ((52l) . the 
correlations between SDMEs are taken into account 
(see Tabs.[15llll). 

For coherent production on the deuteron (iso- 
spin zero), only isoscalar meson exchange is allowed; 
hence there are no contributions from tt, ai, or 61 ex- 
change. 

The incoherent contribution to the cross section on 
the neutron is expected to have an unnatural-parity- 
exchange contribution similar to that for the proton. 
The resulting value of ui for the deuteron is hence 
expected to be smaller than that for the proton due 
to the admixture of coherent scattering. A possible 
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16. The and t' dependence of -ui = 1 
— 2rJi — 2ri_i for proton (filled squares) and deuteron 
(filled circles) data. The values of Mi for yields integrated 
over the range 1 GeV'^ < < 7 GeV^ are shown as 
open symbols. The inner (outer) error bars represent the 
statistical (total) uncertainties. 



indication of this behavior is observed in the lowest 
t' bin of the right section of Fig. [TBI where ui of the 
proton exceeds that of the deuteron. 

The HERMES results on ui,U2, and are pre- 
sented in Fig. [17] and in the top section of Tab. [11] 
The value of is measured here for the first time. The 
combination of proton and deuteron data shows the ex- 
istence of UPE amplitudes on the level of almost three 
standard deviations in the total uncertainty: Ui'^'^ = 
0.106 ± 0.036. In addition, results on ui and U2 from 
other experiments are given in Fig. [17] and in the bot- 
tom panel of Tab.[TTl While U2 is measured to be com- 
patible with zero by all experiments, ui is found to be 
consistent with zero only for high values of VF, as ex- 
pected for TT, ai, and bi exchanges. For low values of 
W, the averaged result from the older measurements, 
yiowW _ g yg _|_ Q ig agreement with the conclu- 
sion that UPE amplitudes exist at HERMES kinemat- 
ics. 

It is worth recalling that the existence of unnatural- 
parity exchange in p'^production by a virtual photon, 
with longitudinally polarized beam and target, results 
in a double-spin asymmetry [21]. At HERMES [l^ this 
asymmetry was found to be non-zero for the proton, 
with a significance of 1.7 standard deviations of the 
total uncertainty; the asymmetry for the deuteron was 
smaller, as discussed in Refs. [T4 1 I2T ] . 

We note that there is no agreement between the 
HERMES measured value of ui at = 3 GcV^ and 
values of ui calculated in variants of a GPD-based 
model 1551. 



10 Contribution of the Helicity-Flip and 
UPE Amplitudes to the Full Cross Section 

Non-conservation of s-channel helicity arises from the 
existence of non-zero helicity-single-flip and/or helicity- 
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Fig. 17. Average values of iti, U2 and 113 calculated ac- 
cording to (|52I54|) from HERMES proton (filled squares) 
and deuteron (filled circles) SDMEs are shown together 
with the values calculated from published SDMEs from 
DESY EZ], SLAC'79 [55], SLAC74 [56], ZEUS BPS [9], 
ZEUS DIS [To] and HI [H]- For HERMES (other exper- 
iments) systematic uncertainties are combined in quadra- 
ture with (without) accounting for correlations between the 
SDMEs. The HERMES deuteron and SLAC'74 data points 
are presented with a small horizontal offset to improve their 
visibility. 



double-flip amplitudes. It can be quantified by mea- 
suring ratios Tij , of helicity- flip amplitudes T^- to the 
square root of lihe sum of all amplitudes squared. 



with TV = eWi 
squared ratio 



(58) 



+ Nt as defined in Section 3. The 
represents the fractional contribution 
from the amplitude T^j to the full cross section. The 
Tij's can be expressed in terms of SDMEs, as shown in 
Appendix C. 

For the helicity- flip amplitude Tqi, describing the 
transition 7^ — *■ p^, the quantity tqi is approximated 
as: 



TOl 



(59) 



For the helicity- flip amplitude Tio, describing the tran- 
sition 72 p^, the quantity tio is given by 



no 



y^(rfi + Im{rf_ J)2 + (Im{r7_ J - rfj^ 



, (60) 
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Fig. 18. Ratios of certain helicity-flip amplitudes to the 
square root of the sum of all amplitudes squared: roi, no, 
and ri_i. HERMES results on proton (filled squares) and 
deuteron (filled circles) are calculated according to (|59ll6ip . 
while results from DESY [33, SLAC'79 [55], SLAC'74 [56], 
ZEUS BPS 0, ZEUS DIS [10] and HI [TT] are calculated 
according to (|62)I64|) . For HERMES (other experiments) 
systematic uncertainties are combined in quadrature with 
(without) accounting for correlations between the SDMEs. 
The HERMES deuteron and SLAC'74 data points are pre- 
sented with a small horizontal offset to improve their visi- 
bility. 



and for the hclicity-doublc-flip amplitude T] 
scribing the transition 7* 
given by 



1-1: 



pj^, the quantity ti 



de- 
_i is 



Tl-l 



,J{rhr + {lni{rU}r 



1-1 



Im{rf_J 



(61) 



The resulting values for proton and deuteron data 
are presented in Fig. [18] and in the top section of 
Tab. [T^ for the entire HERMES kinematic region. Non- 
conservation of s-channel helicity is clearly observed 
for the amplitude Tqi and, for the first time, for the 
amplitude Tio, although with somewhat lower statis- 
tical significance. 

Polarized SDMEs cannot be determined from col- 
lider data, as the collider kinematic conditions imply 
e w 1. According to ([39]) . this suppresses the contribu- 
tion of polarized SDMEs to W^. 



In Ref. [5] the amplitude ratios are approximated 
as follows: 



no 



Tl-l 



' 00 



^04 ' 



V2^oo 
Re{r?4} + Mr!o} 



' 111 



(62) 
(63) 
(64) 



In contrast to (|59ll6ip , these expressions rely on the as- 
sumption of zero phase difference between the consid- 
ered amplitude (Tqi, Tiq, or Ti_i) and the correspond- 
ing dominant amplitude (Too or Tn). Results for the 
quantities Ty from ZEUS and other experiments, cal- 
culated from unpolarized proton SDMEs, are shown in 
Fig. [T5]and in the bottom section of Tab. [T^] 

The combined effect of s-channcl helicity non-con- 
servation and of a contribution of UPE to the full cross 
section can be estimated according to P2I331 127I28P as 
follows. First note that 



da full 



dt 



mil 



(65) 



where 



|Tn|2 (66) 

contains only the contributions of s-channel helicity 
conserving NPE amplitudes. The s-channel helicity 
non-conserving fractional contribution of NPE amp- 
litudes to the cross section is defined as 



(2e|r, 



101 



|Toi 



'01 



\Ti-i\')/Af, 

(67) 

The HERMES result for t| is 0.025 ± O.OOS^tat ± 
Om'isyst and 0.028 ± 0.002^tat ± 0.002^2^,4 for the pro- 
ton and deuteron, respectively. 

Correspondingly, the UPE contribution is defined 

as: 

'^UPE = 

5^(2e|C/ioP + iC/oiP + iC/i-iP + iC/iiH/AAo. (68) 

Because the contributions of amplitudes Uoi and C/i_i 
to ([55]) are negligibly small, T^p^ and ui ([55]) can be 
approximately related to one another as: T^p^ ~ wi/2. 
Accordingly, the first determination of the fractional 
UPE contribution to the full cross section t, 
0.063 ± O.Oll^tat ± 0.025,y,t 



and 0.046 ± 0.008. 



UPE 



IS 
± 



0.023sj,st for the proton and deuteron, respectively. 



11 Longitudinal-to- Transverse Cross 
Section Ratio 

In principle the longitudinal-to-transverse cross sec- 
tion ratio R p4p can experimentally be determined 
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directly from the two cross sections if they can be ex- 
tracted separately from the data using, e.g., the Rosen- 
bluth decomposition technique [57j . For given values 
of and W (or and Xb), this requires data sets 
at different values of e, so that measurements at dif- 
ferent beam energies are necessary No data on 
vector meson production using such a decomposition 
have been published. 
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11.1 Approximations for R 

A common approximation to the ratio R is experimen- 
tally determined from the measured SDME r\ 



,04. 
00- 



' 00 



e 1 



^04 ■ 



(69) 



The quantity R^'^ represents the ratio of cross sections 
for longitudinal and transverse p° polarization, and it 
is not identical to the true R that represents the ratio 
of the cross sections with respect to the polarization 
of the virtual photon. The relation between R^^ and 
R is obtained by comparing (|69l77p with (|34l28l27p : 



e(l + ^) 



(70) 



with 



(1 + ei?04) 



N 

x^il^oil' + M - 2<|T^iol' + |f/ioP)} (71) 

(see Appendix D). In the case of SCHC, ry = and 
R°^ = R. The quantity R°^ can be either smaller or 
larger than i?, depending on the sign of the small pa- 
rameter 77. The latter can be calculated from data by 
neglecting the small contributions of the helicity-flip 
UPE amplitudes C/10, ?7oi in ((7T|) : 



,7«(l + 6i?"^)(ro^i-2er^). 



(72) 



where tqi and tiq are given in (|59ll60p . 

Regge phenomenology suggests that contributions 
of unnatural-parity exchange are more significant at 
the lower energies typical of this experiment, and de- 
crease at collider energies. In order to allow a com- 
parison of HERMES results on R with those at high 
energy and also with GPD-based calculations, the ra- 
tio R^^^ is determined from by subtracting the 
contributions of all UPE amplitudes. The dependence 
of the difference AR^^^ = R°^ - R^^^ on \Uij\^ can 
be determined in a linear approximation as 
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Fig. 19. dependence of the longitudinal-to-transverse 
cross section ratio measured at HERMES. Results from 
proton (deuteron) data are shown in the left (right) panel. 
Filled symbols represent the value of 7?"* calculated from 
rfjQ (|69p . open symbols correspond to the true value of 
R calculated according to (|70l72p . and crosses (diamonds) 
(1751). Total uncertainties are shown, cal- 



represent R 

culated by combining the statistical and systematic uncer- 
tainties in quadrature. The data points for R and R^^^ 
are presented with a small horizontal offset to improve their 
visibility. 



Assuming the hierarchy ([57]) of UPE amplitudes, this 
can be approximated by retaining only Un: 



AR 



UPE 



d\U: 



(73) 



According to ([70]) and (|7ip . i?*^"* can be approximated 
by i? = Nl/Nt, which yields, along with (|27l28p . 



AR 



UPE 



-R 



Mt + eA^L 

Ui 



(l + ei?"4) 



(74) 



Here ui w 2^|C/ii [^/(A/r + (-Ml) is used instead of 
([55)1 . The final approximate formula for R^^^ = R^'^— 
AR^P^ is 



R 



NPE 



^"'[l + ^(l + e^"')] 



(75) 



11.2 HERIVIES Results on R 



Evaluations of R from HERMES data are performed 
for the enthe interval GeV^ < -t' < 0.4 GeV^ 
The dependences of the quantities R'^^ ((69l) and 
R (|70|72p are presented in Fig. [HI In the HERMES 
kinematic conditions, at e « 0.8, the value of 77 is about 
0.1 (—0.1) for the proton (deuteron), and the magni- 
tude of the difference between R and R^^ is small, of 
the order of 0.1. 



25 



□ CLAS <W)=2.1 GeV 

□ CLAS <W)=2.6 GeV 

o CORNELL <W)=3.5 GeV 
■ HERMES <W)=4.8 GeV 
* E665 <W)=18GeV 

A ZEUS <W)=75 GeV 
^ H1 <W)=75 GeV 




Pi 



& 
a 
u 
H 



DJ3 
S 
O 



■ R^™ HERMES 
A R^"* ZFTTS 

-FitofR^™, R«^ 
Ro GPDs 




... 1 




1 


(GeV^) 



(GeVO 

Fig. 20. dependence of the longitudinal-to-transverse cross section ratio for exclusive production on the proton. 
Left panel: i?"* calculated from the SDME roo according to (|69p . HERMES proton data (filled squares) are compared to 
measurements of CLAS [551155] . Cornell [Sg, E665 [BT], HI [TT], and ZEUS PfTD] . The more recent CLAS data [51] (small 
squares) are from a narrow bin m xb with approximately the same {xb) as the HERMES data, which are integrated 
over the xb acceptance. Right panel: B^* for ZEUS (triangles) and R^^^ for HERMES (squares), fitted separately 
according to (|76p . For all data points, total uncertainties are shown. Theoretical calculations [35] of Ro = iTbol^/lTii ^ 
are shown as dashed line at W= 5 GeV; the uncertainties arising from the uncertainties in the parton distribution 
functions are shown as a shaded band 1351. 



In section [10] it was shown that by analyzing the 
amphtudes that comprise the SDMEs, a statistically 
significant, non-zero UPE contribution to the cross 
section exists. At the intermediate energy of the HER- 
MES experiment this contribution is small. If it is 
caused by exchange of tt, ai, or 6i, this contribution 
would be negligible at higher energies [S]. In order to 
compare the HERMES results on R with those of ex- 
periments at higher energy, it is appropriate to correct 
i?""* for the UPE contribution and consider i?^^^. The 
value of AR^P^ is about -0.11 (-0.08) for the pro- 
ton (dcuteron) at HERMES kinematics. The resulting 
values of R^^^ are shown in Fig. [TOland Tab. [T3] 



11.3 Comparison to World Data and Models 

Results for R from different experiments can be com- 
pared only if either R is independent of t', or the t' de- 
pendences of the cross sections and and the t' 
intervals of the measurements of R arc the same. The 
t' dependence of R is determined essentially by the t' 
dependence of the SDME r[j^ (see ([77])), which is found 
to be approximately flat in t' both at HERMES (sec 
Fig. [TO]) and at HI [TT] and ZEUS [TOI kinematics. For 
this case, the ratio of the total cross sections coincides 
with the ratio of the cross sections that are differential 
in t (see (|M|). 



The left panel of Fig. [20] shows HERMES results 
on the dependence of i?""', as measured on the pro- 
ton, in comparison to world data. Given the experi- 
mental uncertainties, there is no discrepancy with the 
data at lower energies from CLAS [5M5^ and COR- 
NELL [60]. The HERMES data at intermediate en- 
ergies are not expected to agree exactly with those 
at high energies because of the UPE contributions ob- 
served in the HERMES data, as discussed in sections[9] 
and [TO] We note that SCHC violating amplitudes are 
also observed in the new CLAS data [SS]. Additional 
reasons may be the importance of valence-quark ex- 
change for NPE amplitudes and also a generally differ- 
ent W dependence of the longitudinal and transverse 
cross sections, as recently discussed in Ref. [35] in the 
context of a GPD-based model. 

The right panel of Fig. [20] presents the HERMES 
results on the longitudinal-to-transverse cross section 
ratio R^^^ , which is corrected for the UPE contribu- 
tions shown in the previous section to be of substantial 
size at intermediate energy. The HERMES data are 
compared to the recent high energy data on i?"^ from 
ZEUS [TO] , for which the UPE contribution is expected 
to be strongly suppressed. 

In order to investigate a possible W dependence of 
the longitudinal-to-transverse cross section ratio, the 
HERMES and ZEUS data are fitted separately to a 
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Q2 dependence suggested by VMD models Pl[71[5^: 

RiQ^) = CO (1^) ' , (76) 

where cq and ci are free parameters and My is the 
mass of the p° meson. The fit resuhs are co = 0.56 ± 
0.08, ci = 0.47 ± 0.12 for HERMES and cq = 0.69 ± 
0.22, ci = 0.59 ±0.15 for ZEUS, with x^/d.o.f. = 0.45 
and 0.15 respectively. These values indicate that 
the fits are dominated by systematic uncertainties. 

A W dependence of the slope is consistent with 
recent calculations using a GPD-based model [35] . We 
note the agreement of these calculations performed at 
W = 5 GeV for values down to 3 GeV^ (see dashed 
curve in Fig. [20]) with the highest — 3 GeV^ point 
of HERMES. Uncertainties in the model calculations 
originating from uncertainties in the parton distribu- 
tions employed are shown as a shaded band superim- 
posed on the curve. 



12 Summary 

HERMES has studied exclusive p° production on the 
proton and deuteron at intermediate energies {{W) = 
4.8 GcV), at the average values of (Q^) = 1.95 GeV^, 
(-f) = 0.13 GeV^, and {xb) = 0.08, using polarized 
beams and unpolarizcd targets. 

By performing a maximum likelihood fit, fifteen 
unpolarizcd SDMEs and, for the first time, eight po- 
larized SDMEs are obtained. The measured SDMEs 
are grouped according to their theoretically expected 
hierarchy. This facilitates the investigation of the rela- 
tive importance of various hclicity amplitudes describ- 
ing different 7* — > transitions. Within the given 
experimental uncertainties, the expected hierarchy of 
relative sizes of hclicity amplitudes is observed. 

Non-zero values are observed for the two helicity- 
flip amplitudes Tqi and Tio, indicating a small but 
statistically significant deviation from the hypothesis 
of s-channcl hclicity conservation. 

The phase difference between the helicity-conser- 
ving amplitudes Tn and Tqo is confirmed to be signif- 
icantly non-zero and is also seen to have a possible 
dependence. For the first time, the sign of the phase 
difference is determined using the polarized SDMEs. 

The kinematic dependences of all 23 SDMEs are 
measured for both hydrogen and deuterium targets. 
Clear dependences on and t' are observed for cer- 
tain SDMEs. No significant difference between proton 
and deuteron results is seen. 

The evaluation of certain linear combinations of 
SDMEs provides an indication that at the interme- 
diate energy of the HERMES experiment, contribu- 
tions of unnatural-parity-exchange amplitudes exist. 
Such amplitudes are naturally generated by a quark- 
exchange mechanism. 



In order to determine the longitudinal-to- 
transverse cross section ratio with respect to the 
polarization of the virtual photon, an approximation 
to the ratio R of the cross sections for longitudinal 
and transverse p° polarizations is calculated from the 
SDME function of Q^. The results obtained 

for other SDMEs permit us to improve this approx- 
imation of R by taking into account transitions of 
natural parity that do not conserve s-channel helicity. 
In order to facilitate comparison with high-energy 
collider data, a correction is applied to to exclude 
contributions from unnatural-parity exchange. The 
comparison of the dependences of R at low and 
high values of W suggests a possible W dependence 
of the ratio. 
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Appendix A. The 23 Spin-Density Matrix 
Elements Expressed in Terms of Helicity 
Amplitudes 

The basic expressions for the 23 spin density matrix el- 
ements measurable with a polarized lepton beam and 
an unpolarizcd target, ordered according to the ex- 
pected hierarchy of amplitudes, are given in ([77 |) -([99 |) . 
(COMMENT FOR THE EDITORS: these formulae 
appear on the next page). 

Appendix B. Derivation of Formulae for cos S 

and sin 6 

Neglecting small contributions of the products of spin- 
flip amplitudes TwT*^ and UiqU^^ in dSDll-dMl) we ob- 
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A : 7I Pl and -y^ ^ pr 

rot = J2i^\no\' + \Toi\' + \Uoi\'/}/^, (77) 

rl-i = lj2i\Tiif + iTi-il' - \Uiif - \Ui-if}/Af, (78) 

Im{r?_i} = + |Ti-i|" + \Uiif - |(7i-i|'}/AA, (79) 

B : interference of 7^ — * pi and 7^ ^ Pt 

Re{rlo} = -^^Re{2rioroi + (Tn - Ti_i)T,ro}/AA , (80) 

Im{r?o} = ^^Re{2C/io[/oi - (Tn + Ti_i)Too}/AA , (81) 

Im{rIo} = -i=^Im{2[/io;7o*i + (Tn + Ti_i)Too}/AA , (82) 

Re{r?o} = -^^Im{-2rioroi + (Tn - Tj-j)Too}/Af , (83) 

C* 
: iT^ Pl 

M^?o} = ^ReieTioToo + iroi(rii - Ti-i)* + ic/oi(;7ii + )7i-i)*}/A/' , (84) 

ReKo} = i^Re{-Toi(rii ~ Ti-i)* + C/oi(f/ii + f/i-i)*}/AA, (85) 

Im{r?o} = i^Re{roi(rii +ri-i)* - Um{Uii - Ui-i)*}/M, (86) 

rj^o = \/2^Re{ToiT,ro}/AA , (87) 

Im{r?o} = -i^Im{roi(rii +ri_i)* + ?7oi(;7ii - f/i-i)*}/AA, (89) 

r-oo = V2^Ini{roiroo}/AA , (90) 
D : iL-^ Pt 

r\i = --^^Re{Tio(Tii - Ti-i)* + ?7io(;7ii - Ui-^)*}/U, (91) 

r5^-i = ^^Re{-rio(rii - Ti-i)* + ;7io((7ii - t/i-i)*}/AA, (92) 
V 2 

Im{r?_i} = --L^Re{rio(Tii + Ti-i)* - t/io(f/ii + ;7i_i)*}/AA, (93) 

Im{rLi} = ^^Im{Tio(Tn +Ti_i)* - Uw{Uii + ;7i-i)*}/AA, (94) 
V 2 

r?i = — ^^Im{Tio(Tii - Ti-i)* + t/io(;7ii - ;7i_i)*}/AA, (95) 
v2 

r?_i = -^^Im{rio(rii -Ti-i)* -f/io(i7ii -f/i-i)*}/AA, (96) 

E : 7* T — * Pt 

= ^Re{-6!rio|" + e\Uw? + Ti-iTn - U^-^m^}/U, (97) 

rli = ^Re{Ti_iri*i + U^-^Ul,}/N , (98) 

Im{r?_i} = -^MTi-iTii - Ui-im,}/Af , (99) 
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tain the approximate relations: 

V2{Re{rlo} - li-n{r%}) « Re{Tnn,)/M 

= iTiillToolcosVAA, (100) 
V2{Re{r%} + Im{r[„}) « Im(TnT*„)/AA 

^\Tn\\Too\smS/Af . (101) 

Neglecting |Toip and |i7oiP in the numerator in rela- 
tion (|77)) we get 



(102) 



Reealling the formulae (gT]), (US]), and (gH]) for TV, we 
obtain from (j77p the approximate relation 



' 00 



|Tiip + |[/np /A^ (103) 



if we neglect small contributions of spin-flip amplitudes 
in the numerator of (|103[) . As seen from relation (|103p 
the difference 1 — Tqq contains |C/iiP; to cancel this 
contribution, the difference of two SDMEs defined by 
(l78l) and (17911 is considered: 



rl_,-lm{rU}=Y,[\Tr 



\Uu 



lAf . (104) 



Adding (|T03l) and ([T04l) together the relation 



(1 



Im{r?_i})/2«|rn|VA^ (105) 



is obtained. Dividing pOOp and (jlOip by the square 
root of the product of (|102p and (|105p . formulae ([ST 
and ([5T|) are obtained respectively. 



Appendix C. Derivation of Formula for tio 

Here we derive the formula for tio only; formu- 
lae for tqi and ri_i can be derived in an analo- 
gous way. If we retain only linear contributions of 
small s-channel helicity-ffip amplitudes in the basic 
formulae for SDMEs jfTj-jMl), neglecting unnatural- 
parity-exchange contributions and bilinear products of 
helicity-flip amplitudes we obtain: 



'00 



elToolVA/'o , 



(106) 



Rc{r?on = Rc{eTioT*o + -T^^T;^]/No , (107) 

M'^i'o} = -^Rc{roir*i}/AAo , (i08) 

ri_i = -Im{rti} = \\Tii?/Mo , (109) 

= Im{rti} = i=Re{TioT*i}/AAo , (110) 



Im{r[_i} 



' 11 



V2' 
1 
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Im{rioTi*J/AAo , (111) 



where A/n = e|T} 



00 I 



iTiip (see dSni)). 



The parameter tiq was defined in [9] by the relation 

no = |Tio|/v/|rooP + |TnP (112) 
and estimated from (|106ll07ll08p with the formula 
Re{r?cJ}+Re{rio} 



Tio 



V ' 00 

iTiol cos Jio 



l^iiPA 



(113) 



where 5io = arg(Tio) — arg{ToQ). Comparison of 
pi2p and (jll3p shows that they are equal to each other 
if e « 1 and 5io w 0. 

Instead of (|113p we derive a formula which is ap- 
plicable for any values of 6io and e. Combining (jllOp 
and piip . we obtain 



^ f(Re{TioT*i})2 + (ImlTioTAD^ 



2 
A^ 



Tio|'|Tn|2. 



(114) 



Dividing (UHl) by 2rJ_i - 2Im{r^_J = 2\Tn\^-/Afo 
(see (|109p ) we get the final approximate formula 



no 



\Tio\/VNo 



irl,+lm{rU}r + ilm{rl_,} 



2{r\_, - Im{rti}) 



(115) 



Since Ao = A/" within the approximation under consid- 
eration, formula (|115p corresponds to definition (|58p of 
TiQ. In the case e = 1, which was considered in Ref. [9], 
the estimate (|115l) for tiq coincides with the definition 
given in pi2p . 



Appendix D. Derivation of Relations 
between R"'* and R 



From ()69I77|26P it follows that 

EWT^ooP + lToiP + lC/oiP} 



.04 



Nt + eNL - E{e|rooP + iToiP + \Uoi\^} 
^^L + E{-2e(|rioP + \Ui^?) + iToiP + \U^i?} 

Mt - E{-2e(|rioP + |J7ioP) + iToiP + |i7oiP} ' 

(116) 

where in the second step we have used the formula for 
Al (|28p for the transformation of (|116p . Dividing both 
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the numerator and denominator in (jll6p by A/r and 
remembering the definition p4p of R we get 



ei? + C 



with 



C = ^{-2e(|TioP + iC/ioH + \Toi\^ + |C/oi|'}/A/r 
_ 77(1 + ei?) 



(118) 



The last relation follows from comparison of pisp and 
the definition (|7T|) of rj. Equation (|117p can be easily 
rewritten in the form 



e 



which is equivalent to ((70|) if we take into account the 
relation pisp between 77 and (. 



Appendix E. Kinematic Intervals, Mean Val- 
ues for Kinematic Variables and SDMEs, for 
Proton and Deuteron 

The resulting SDMEs with statistical and systematic 
uncertainties are presented below in tabular form for 
hydrogen and deuterium targets. First, in Tab. [2] the 
mean kinematic values are presented for the entire 
kinematic region of the measurement and for each bin 
used in the Q^, t' and xb dependences. In Tabs. [3l 
HI E] (in [ll [8|) the results of the measurement of the 
Q^, t' and xs-dependences, respectively, for the pro- 
ton (deuteron) are listed. The values of the phase dif- 
ference S between Tn and Tqq amplitudes, from proton 
and deuteron data, are contained in Tab. [51 The kine- 
matic dependences of the ui value used for the study of 
unnatural-parity-exchange amplitudes are in Tab. 1101 
The SDMEs measured over the entire kinematic re- 
gion from proton data, but presented using the recent 
representation of Ref. [25] are listed in Tab. [TH The 
correlation matrices of the 23 SDMEs measured for the 
proton and deuteron over the entire kinematic region 
are presented in Tabs. [T5l and [TBI 
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Table 2. The definition of intervals and tiie mean values for kinematic variables for hydrogen (deuterium) data. 



bin 




1 GcV^ < Q2 < 7 GeV^ 


1.95 (1.94) 


0.5 GeV^ <Q^ < 1.0 GeV^ 
1.0 GeV^ < < 1.4 GeV^ 
1.4 GeV^ < Q2 < 2 GeV^ 
2 GeV^ < < 7 GeV^ 


0.82 (0.82) 
1.19 (1.18) 
1.66 (1.66) 
3.06 (3.04) 




{t'), GeV" 


0.0 GeV^ < -t' < 0.04 GeV^ 
0.04 GeV^ < -t' < 0.10 GeV^ 
0.10 GeV^ < -t' < 0.20 GeV^ 
0.20 GeV^ < -t' < 0.40 GoV^ 


0.019 (0.018) 
0.068 (0.068) 
0.146 (0.145) 
0.281 (0.283) 




(x-s) 


0.0 <XB < 0.05 
0.05 <XB < 0.08 
0.08 <XB < 0.35 


0.042 (0.042) 
0.064 (0.064) 
0.120 (0.119) 


1 GcV^ <Q'' <7 GeV^ 


(e) = 0.80 ±0.01 



Table 3. The 23 unpolarized and polarized SDMEs for p" production from the proton in bins defined by the limits 
0.5, 1.0, 1.4, 2.0, and 7.0 GeV^ The first uncertainties are statistical, the second systematic. 



Element 


(Q^) = 0.82 GeV^ 


(Q^) = 1.19 GeV^ 


(Q^) = 1.66 GeV^ 


= 3.06 GeV^ 


r°ot 




0.349 ± 0.026 ± 0.061 


0.368 ± 0.018 ± 0.011 


0.397 ± 0.017 ± 0.018 


0.454 ± 0.014 ± 0.011 


rl- 




0.283 ± 0.023 ± 0.049 


0.262 ± 0.018 ± 0.024 


0.274 ± 0.019 ± 0.024 


0.204 ± 0.017 ± 0.012 


Im 




-0.294 ± 0.019 ± 0.038 


-0.255 ± 0.016 ± 0.022 


-0.239 ± 0.017 ± 0.011 


-0.197 ± 0.017 ± 0.012 


Re 


rlo 


0.151 ± 0.028 ± 0.026 


0.171 ± 0.007 ± 0.000 


0.161 ± 0.006 ± 0.004 


0.141 ± 0.006 ± 0.008 


Im 


rfo 


-0.149 ± 0.015 ± 0.010 


-0.167 ± 0.007 ± 0.003 


-0.167 ± 0.006 ± 0.005 


-0.156 ± 0.006 ± 0.010 


Im 


rlo 


0.079 ± 0.068 ± 0.011 


0.092 ± 0.038 ± 0.010 


0.039 ± 0.036 ± 0.004 


0.187 ± 0.034 ± 0.018 


Re 


rfo 


0.040 ± 0.043 ± 0.011 


0.020 ± 0.031 ± 0.008 


0.074 ± 0.034 ± 0.002 


0.098 ± 0.032 ± 0.005 


Re 


r^o 


0.028 ± 0.028 ± 0.020 


0.029 ± 0.007 ± 0.003 


0.035 ± 0.007 ± 0.011 


0.026 ± 0.007 ± 0.003 


Re 


rlo 


-0.037 ± 0.044 ± 0.032 


-0.043 ± 0.012 ± 0.006 


-0.036 ± 0.012 ± 0.012 


-0.009 ± 0.013 ± 0.010 


Im 


rlo 


0.023 ± 0.019 ± 0.007 


0.022 ± 0.012 ± 0.018 


0.005 ± 0.012 ± 0.024 


0.022 ± 0.013 ± 0.008 


''00 




0.121 ± 0.038 ± 0.039 


0.094 ± 0.017 ± 0.017 


0.057 ± 0.015 ± 0.019 


0.151 ± 0.015 ± 0.007 


roo 




-0.054 ± 0.039 ± 0.013 


0.011 ± 0.032 ± 0.018 


0.007 ± 0.031 ± 0.009 


0.037 ± 0.034 ± 0.002 


Im 


rlo 


0.002 ± 0.041 ± 0.008 


-0.041 ± 0.026 ± 0.005 


-0.074 ± 0.025 ± 0.005 


0.048 ± 0.024 ± 0.006 


''oo 




0.022 ± 0.079 ± 0.026 


0.040 ± 0.084 ± 0.014 


0.054 ± 0.086 ± 0.011 


0.010 ± 0.085 ± 0.016 


rl^ 




-0.015 ± 0.010 ± 0.007 


-0.011 ± 0.006 ± 0.006 


-0.008 ± 0.006 ± 0.011 


-0.021 ± 0.006 ± 0.016 


rl- 


1 


0.009 ± 0.011 ± 0.019 


0.008 ± 0.007 ± 0.006 


-0.013 ± 0.007 ± 0.003 


0.020 ± 0.007 ± 0.008 


Im 




-0.011 ± 0.010 ± 0.013 


0.002 ± 0.007 ± 0.007 


0.002 ± 0.007 ± 0.004 


-0.010 ± 0.007 ± 0.007 


Im 


rU 


-0.003 ± 0.078 ± 0.021 


0.023 ± 0.056 ± 0.013 


-0.005 ± 0.055 ± 0.010 


-0.109 ± 0.047 ± 0.004 


r?i 




0.019 ± 0.053 ± 0.007 


0.056 ± 0.045 ± 0.004 


0.051 ± 0.044 ± 0.006 


-0.002 ± 0.035 ± 0.005 


rl^ 


1 


0.013 ± 0.062 ± 0.008 


0.072 ± 0.053 ± 0.011 


-0.018 ± 0.054 ± 0.004 


0.004 ± 0.045 ± 0.014 


r?i 


1 


-0.024 ± 0.013 ± 0.021 


-0.014 ± 0.010 ± 0.010 


-0.019 ± 0.010 ± 0.003 


0.001 ± 0.009 ± 0.007 


r\. 




-0.039 ± 0.017 ± 0.018 


-0.034 ± 0.013 ± 0.013 


-0.023 ± 0.013 ± 0.008 


-0.018 ± 0.012 ± 0.010 


Im 


rU 


0.021 ± 0.051 ± 0.010 


0.000 ± 0.033 ± 0.004 


-0.031 ± 0.032 ± 0.007 


-0.026 ± 0.028 ± 0.005 
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Table 4. The 23 unpolarized and polarized SDMEs for production from the proton in —t' bins defined by the hmits 
0.0, 0.04, 0.1, 0.2, and 0.4 GeV^. The first uncertainties are statistical, the second systematic. 



Element 


{-t') = 0.019 GeV^ 


(-<') = 0.068 GeV^ 


= 0.146 GeV^ 


{-t') = 0.281 GeV ^ 






0.393 ± 0.018 ± 0.019 


0.394 ± 0.018 ± 0.025 


0.415 ± 0.019 ± 0.021 


0.481 ± 0.019 ± 0.028 


r\- 


1 


0.235 ± 0.021 ± 0.014 


0.265 ± 0.021 ± 0.012 


0.260 ± 0.020 ± 0.031 


0.203 ± 0.020 ± 0.030 


Im 


rl-i 


-0.204 ± 0.020 ± 0.008 


-0.243 ± 0.019 ± 0.011 


-0.232 ± 0.019 ± 0.041 


-0.231 ± 0.019 ± 0.021 


Re 


rlo 


0.156 ± 0.007 ± 0.005 


0.156 ± 0.007 ± 0.010 


0.153 ± 0.007 ± 0.007 


0.153 ± 0.008 ± 0.009 


Im 


r% 


-0.162 ± 0.007 ± 0.005 


-0.170 ± 0.007 ± 0.010 


-0.155 ± 0.007 ± 0.006 


-0.153 ± 0.008 ± 0.008 


Im 


rlo 


0.103 ± 0.040 ± 0.007 


0.112 ± 0.040 ± 0.007 


0.081 ± 0.042 ± 0.011 


0.163 ± 0.047 ± 0.033 


Re 


rlo 


0.042 ± 0.037 ± 0.011 


0.059 ± 0.037 ± 0.006 


0.100 ± 0.036 ± 0.013 


0.114 ± 0.039 ± 0.015 


Re 


rlt 


0.018 ± 0.008 ± 0.004 


0.027 ± 0.008 ± 0.011 


0.035 ± 0.008 ± 0.007 


0.038 ± 0.008 ± 0.003 


Re 


r\o 


-0.009 ± 0.014 ± 0.005 


-0.045 ± 0.014 ± 0.024 


-0.013 ± 0.015 ± 0.009 


-0.046 ± 0.016 ± 0.015 


Im 


rlo 


-0.001 ± 0.014 ± 0.011 


0.030 ± 0.014 ± 0.019 


0.015 ± 0.013 ± 0.019 


0.030 ± 0.015 ± 0.013 


rlo 




0.039 ± 0.016 ± 0.001 


0.068 ± 0.016 ± 0.030 


0.136 ± 0.018 ± 0.010 


0.219 ± 0.020 ± 0.022 


'"oo 




0.019 ± 0.035 ± 0.020 


0.015 ± 0.035 ± 0.027 


-0.026 ± 0.036 ± 0.021 


-0.005 ± 0.041 ± 0.013 


Im 


rlo 


-0.035 ± 0.028 ± 0.001 


-0.044 ± 0.027 ± 0.004 


0.018 ± 0.029 ± 0.007 


0.009 ± 0.032 ± 0.025 


''00 




-0.013 ± 0.103 ± 0.011 


0.128 ± 0.097 ± 0.012 


-0.028 ± 0.095 ± 0.028 


0.066 ± 0.097 ± 0.012 






r\ r\r\r\ i r\ r\r\^ i r\ r\r\ a 

—0.009 ± 0.007 ± 0.004 


r\ r\'i 1 1 r\ r\i~\f^ i r\ r\r\r\ 

—0.011 ± 0.006 ± 0.009 


r\ r\r\r\ i r\ r\r\^ i r\ r\'\ t 

—0.009 ± 0.007 ± 0.011 


—0.031 ± 0.006 ± 0.024 


A- 


1 


-0.009 ± 0.008 ± 0.003 


0.005 ± 0.008 ± 0.004 


0.003 ± 0.008 ± 0.008 


0.012 ± 0.008 ± 0.008 


Im 




0.010 ± 0.008 ± 0.002 


-0.010 ± 0.008 ± 0.008 


0.006 ± 0.008 ± 0.003 


-0.009 ± 0.008 ± 0.002 


Im 




-0.034 ± 0.064 ± 0.011 


-0.040 ± 0.060 ± 0.004 


-0.076 ± 0.060 ± 0.015 


-0.005 ± 0.058 ± 0.008 


r-fi 




0.018 ± 0.050 ± 0.003 


-0.010 ± 0.048 ± 0.009 


0.061 ± 0.045 ± 0.006 


0.068 ± 0.042 ± 0.006 


rf- 


1 


0.021 ± 0.062 ± 0.005 


0.024 ± 0.058 ± 0.009 


-0.019 ± 0.055 ± 0.007 


0.051 ± 0.053 ± 0.006 


tT- 


1 


0.008 ± 0.011 ± 0.003 


0.008 ± 0.011 ± 0.010 


-0.033 ± 0.011 ± 0.008 


-0.029 ± 0.010 ± 0.003 


rli 




-0.022 ± 0.015 ± 0.016 


0.002 ± 0.015 ± 0.009 


-0.036 ± 0.014 ± 0.003 


-0.034 ± 0.014 ± 0.012 


Im 


rU 


-0.038 ± 0.036 ± 0.008 


-0.015 ± 0.035 ± 0.001 


-0.014 ± 0.036 ± 0.006 


-0.036 ± 0.036 ± 0.010 
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Table 5. The 23 unpolarized and polarized SDMEs for p production from the proton in xb bins defined by the hmits 
0.0, 0.05, 0.08, and 0.35. The first uncertainties are statistical, the second systematic. 



Element 


(xb) =0.042 


{a;fl)=0.064 


{a;s)=0.120 


rlt 




0.349 ± 0.044 ± 0.031 


0.405 ± 0.017 ± 0.013 


0.448 ± 0.014 ± 0.012 


r\- 


1 


0.297 ± 0.036 ± 0.006 


0.253 ± 0.014 ± 0.023 


0.216 ± 0.015 ± 0.019 


Im 


rU 


-0.308 ± 0.032 ± 0.024 


-0.229 ± 0.031 ± 0.008 


-0.198 ± 0.014 ± 0.015 


Re 


rfo 


0.151 ± 0.027 ± 0.014 


0.176 ± 0.008 ± 0.007 


0.146 ± 0.006 ± 0.012 


Im 


r% 


-0.172 ± 0.015 ± 0.006 


-0.169 ± 0.005 ± 0.003 


-0.157 ± 0.006 ± 0.012 


Im 


rlo 


0.132 ± 0.067 ± 0.015 


0.045 ± 0.061 ± 0.016 


0.160 ± 0.032 ± 0.018 


Re 


rfo 


0.021 ± 0.041 ± 0.002 


0.067 ± 0.042 ± 0.000 


0.092 ± 0.032 ± 0.011 


Re 


rll 


0.023 ± 0.021 ± 0.005 


0.038 ± 0.017 ± 0.010 


0.033 ± 0.006 ± 0.005 


Re 


rlo 


-0.013 ± 0.035 ± 0.014 


-0.050 ± 0.011 ± 0.010 


-0.021 ± 0.011 ± 0.019 


Im 


rlo 


0.035 ± 0.021 ± 0.017 


0.009 ± 0.022 ± 0.015 


0.027 ± 0.011 ± 0.017 


roo 




0.079 ± 0.051 ± 0.034 


0.056 ± 0.013 ± 0.017 


0.130 ± 0.013 ± 0.014 


■Too 




0.044 ± 0.064 ± 0.016 


0.000 ± 0.029 ± 0.005 


0.067 ± 0.030 ± 0.007 


Im 


rfo 


-0.025 ± 0.042 ± 0.001 


-0.067 ± 0.033 ± 0.001 


0.014 ± 0.022 ± 0.005 


rlo 




0.100 ± 0.092 ± 0.011 


0.030 ± 0.074 ± 0.003 


-0.005 ± 0.090 ± 0.031 


rli 




0.005 ± 0.022 ± 0.011 


-0.014 ± 0.011 ± 0.005 


-0.019 ± 0.004 ± 0.016 


rl- 




-0.014 ± 0.020 ± 0.002 


0.017 ± 0.007 ± 0.004 


0.000 ± 0.006 ± 0.005 


Im 




0.000 ± 0.017 ± 0.003 


0.003 ± 0.014 ± 0.007 


-0.001 ± 0.006 ± 0.007 


Im 


ri-1 


0.019 ± 0.099 ± 0.012 


0.071 ± 0.135 ± 0.019 


-0.142 ± 0.044 ± 0.010 


rfi 




0.042 ± 0.060 ± 0.007 


0.081 ± 0.079 ± 0.008 


-0.002 ± 0.033 ± 0.009 


rl- 


1 


0.089 ± 0.065 ± 0.003 


0.031 ± 0.077 ± 0.003 


0.005 ± 0.043 ± 0.009 


r?i 


1 


-0.025 ± 0.020 ± 0.007 


-0.001 ± 0.012 ± 0.010 


-0.016 ± 0.008 ± 0.003 


ri\ 




-0.067 ± 0.030 ± 0.008 


-0.018 ± 0.014 ± 0.011 


-0.032 ± 0.010 ± 0.009 


Im 


rU 


-0.020 ± 0.067 ± 0.014 


-0.017 ± 0.061 ± 0.009 


-0.001 ± 0.026 ± 0.004 
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Table 6. The 23 unpolarized and polarized SDMEs for p° production from the deuteron in bins defined by the 
limits 0.5, 1.0, 1.4, 2.0, and 7.0 GeV^. The first uncertainties are statistical, the second systematic. 



SDME 


(Q2) = 0.82 GeV^ 


{Q2) = 1.18 GeV^ 


= 1.66 GeV^ 


(Q^) = 3.04 GeV^ 


rlt 




0.365 ± 0.015 ± 0.058 


0.386 ± 0.017 ± 0.010 


0.445 ± 0.012 ± 0.022 


0.407 ± 0.011 ± 0.013 


r\- 


1 


0.294 ± 0.013 ± 0.060 


0.278 ± 0.016 ± 0.022 


0.235 ± 0.013 ± 0.012 


0.216 ± 0.014 ± 0.008 


Im 


rU 


-0.287 ± 0.014 ± 0.047 


-0.267 ± 0.013 ± 0.034 


-0.194 ± 0.014 ± 0.020 


-0.219 ± 0.014 ± 0.017 


Re 


rfo 


0.159 ± 0.008 ± 0.025 


0.173 ± 0.005 ± 0.002 


0.160 ± 0.005 ± 0.006 


0.154 ± 0.005 ± 0.005 


Im 


r% 


-0.159 ± 0.008 ± 0.015 


-0.165 ± 0.005 ± 0.001 


-0.155 ± 0.005 ± 0.003 


-0.139 ± 0.005 ± 0.007 


Im 


rlo 


0.070 ± 0.042 ± 0.004 


0.080 ± 0.028 ± 0.011 


0.102 ± 0.030 ± 0.004 


0.125 ± 0.026 ± 0.007 


Re 


rfo 


0.080 ± 0.028 ± 0.012 


0.080 ± 0.023 ± 0.003 


0.132 ± 0.026 ± 0.006 


0.117 ± 0.024 ± 0.012 


Re 




0.021 ± 0.007 ± 0.019 


0.036 ± 0.006 ± 0.004 


0.026 ± 0.006 ± 0.016 


0.023 ± 0.005 ± 0.006 


Re 


rlo 


-0.025 ± 0.011 ± 0.031 


-0.017 ± 0.009 ± 0.002 


-0.014 ± 0.010 ± 0.015 


-0.031 ± 0.010 ± 0.013 


Im 


rlo 


-0.006 ± 0.011 ± 0.017 


0.016 ± 0.009 ± 0.018 


0.015 ± 0.010 ± 0.038 


0.004 ± 0.010 ± 0.018 


rlo 




0.097 ± 0.017 ± 0.035 


0.088 ± 0.013 ± 0.011 


0.113 ± 0.012 ± 0.018 


0.119 ± 0.012 ± 0.006 


"Too 




0.019 ± 0.028 ± 0.014 


-0.018 ± 0.024 ± 0.019 


-0.031 ± 0.025 ± 0.020 


-0.036 ± 0.026 ± 0.009 


Im 


rio 


0.005 ± 0.029 ± 0.004 


0.003 ± 0.019 ± 0.008 


0.024 ± 0.021 ± 0.003 


0.062 ± 0.018 ± 0.003 


rlo 




0.138 ± 0.061 ± 0.010 


0.221 ± 0.066 ± 0.021 


0.058 ± 0.069 ± 0.016 


-0.098 ± 0.063 ± 0.007 


rl. 




-0.009 ± 0.006 ± 0.005 


-0.013 ± 0.004 ± 0.009 


-0.017 ± 0.004 ± 0.013 


-0.027 ± 0.004 ± 0.015 


rl- 




0.008 ± 0.007 ± 0.016 


0.009 ± 0.005 ± 0.002 


0.006 ± 0.006 ± 0.004 


0.021 ± 0.006 ± 0.011 


Im 




-0.007 ± 0.007 ± 0.018 


-0.003 ± 0.005 ± 0.004 


-0.003 ± 0.006 ± 0.006 


-0.013 ± 0.006 ± 0.005 


Im 


rU 


-0.066 ± 0.052 ± 0.008 


-0.040 ± 0.040 ± 0.016 


-0.026 ± 0.044 ± 0.001 


-0.100 ± 0.037 ± 0.013 


r?i 




0.007 ± 0.039 ± 0.003 


-0.011 ± 0.035 ± 0.010 


0.047 ± 0.033 ± 0.008 


0.037 ± 0.028 ± 0.001 


rf- 


1 


-0.015 ± 0.047 ± 0.014 


-0.055 ± 0.040 ± 0.015 


-0.083 ± 0.041 ± 0.003 


-0.072 ± 0.036 ± 0.017 


r?i 


1 


0.000 ± 0.009 ± 0.021 


0.003 ± 0.008 ± 0.012 


-0.006 ± 0.008 ± 0.009 


-0.008 ± 0.007 ± 0.003 






-0.029 ± 0.012 ± 0.019 


-0.002 ± 0.009 ± 0.013 


-0.003 ± 0.010 ± 0.011 


-0.012 ± 0.010 ± 0.007 


Im 


rU 


0.006 ± 0.031 ± 0.010 


-0.017 ± 0.024 ± 0.008 


-0.023 ± 0.026 ± 0.003 


0.029 ± 0.022 ± 0.004 
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Table 7. The 23 unpolarized and polarized SDMEs for p" production from the deuteron in —t' bins defined by the 
limits 0.0, 0.04, 0.1, 0.2, and 0.4 GeV^. The first uncertainties are statistical, the second systematic. 



Element 


{-t') = 0.018 Ge-y^ 


(-f') = 0.068 GeV^ 


{-t') = 0.145 GeV^ 


(-*') = 0.283 GeV^ 






0.440 ± 0.014 ± 0.017 


0.396 ± 0.015 ± 0.024 


0.389 ± 0.015 ± 0.017 


0.434 ± 0.017 ± 0.032 




1 


0.225 ± 0.015 ± 0.014 


0.261 ± 0.017 ± 0.011 


0.272 ± 0.016 ± 0.029 


0.247 ± 0.018 ± 0.029 


Im 




-0.208 ± 0.014 ± 0.015 


-0.258 ± 0.015 ± 0.019 


-0.239 ± 0.015 ± 0.028 


-0.229 ± 0.016 ± 0.026 


Re 


rlo 


0.155 ± 0.005 ± 0.006 


0.153 ± 0.006 ± 0.010 


0.172 ± 0.006 ± 0.009 


0.154 ± 0.007 ± 0.009 


Im 


rlo 


-0.154 ± 0.005 ± 0.003 


-0.150 ± 0.005 ± 0.004 


-0.150 ± 0.006 ± 0.006 


-0.154 ± 0.007 ± 0.010 


Im 


rlo 


0.084 ± 0.029 ± 0.006 


0.122 ± 0.028 ± 0.008 


0.108 ± 0.034 ± 0.007 


0.086 ± 0.039 ± 0.015 


Re 


rlo 


0.161 ± 0.026 ± 0.010 


0.095 ± 0.028 ± 0.007 


0.081 ± 0.027 ± 0.006 


0.099 ± 0.030 ± 0.012 


Re 


rlt 


0.016 ± 0.006 ± 0.005 


0.025 ± 0.006 ± 0.011 


0.034 ± 0.006 ± 0.005 


0.046 ± 0.007 ± 0.005 


Re 


rlo 


-0.010 ± 0.010 ± 0.005 


-0.002 ± 0.011 ± 0.021 


-0.020 ± 0.011 ± 0.010 


-0.047 ± 0.013 ± 0.018 


Im 


rlo 


0.012 ± 0.010 ± 0.010 


0.018 ± 0.010 ± 0.018 


0.005 ± 0.011 ± 0.018 


0.036 ± 0.013 ± 0.017 


'-00 




0.053 ± 0.012 ± 0.003 


0.085 ± 0.013 ± 0.033 


0.108 ± 0.014 ± 0.003 


0.215 ± 0.018 ± 0.032 


''oo 




-0.077 ± 0.026 ± 0.025 


-0.001 ± 0.028 ± 0.034 


-0.053 ± 0.029 ± 0.028 


-0.040 ± 0.035 ± 0.030 


Im 




0.020 ± 0.020 ± 0.003 


0.043 ± 0.020 ± 0.003 


0.022 ± 0.023 ± 0.004 


0.033 ± 0.027 ± 0.006 


'-00 




0.157 ± 0.073 ± 0.013 


0.111 ± 0.076 ± 0.011 


-0.147 ± 0.070 ± 0.013 


0.078 ± 0.079 ± 0.012 






—0.018 ± 0.005 ± 0.007 


—0.014 ± 0.005 ± 0.009 


—0.012 ± 0.005 ± 0.005 


—0.025 ± 0.006 ± 0.023 


rl- 


1 


0.001 ± 0.006 ± 0.003 


0.007 ± 0.006 ± 0.004 


0.006 ± 0.006 ± 0.007 


0.017 ± 0.007 ± 0.008 


Im 




-0.007 ± 0.006 ± 0.002 


-0.005 ± 0.006 ± 0.005 


0.001 ± 0.006 ± 0.009 


-0.003 ± 0.007 ± 0.002 


Im 


'^1-1 


-0.037 ± 0.043 ± 0.006 


-0.041 ± 0.045 ± 0.004 


-0.060 ± 0.046 ± 0.002 


-0.080 ± 0.048 ± 0.008 


rl. 




-0.030 ± 0.036 ± 0.002 


-0.003 ± 0.037 ± 0.004 


0.078 ± 0.035 ± 0.001 


0.056 ± 0.037 ± 0.008 


rf- 


1 


-0.095 ± 0.044 ± 0.010 


-0.055 ± 0.044 ± 0.008 


-0.030 ± 0.044 ± 0.008 


-0.091 ± 0.045 ± 0.006 




1 


0.016 ± 0.008 ± 0.004 


-0.005 ± 0.009 ± 0.006 


-0.021 ± 0.009 ± 0.001 


-0.013 ± 0.009 ± 0.008 


'^11 




0.026 ± 0.011 ± 0.013 


-0.004 ± 0.012 ± 0.015 


-0.036 ± 0.012 ± 0.004 


-0.016 ± 0.012 ± 0.011 


Im 




0.003 ± 0.025 ± 0.002 


-0.036 ± 0.026 ± 0.005 


-0.033 ± 0.028 ± 0.002 


0.052 ± 0.031 ± 0.005 
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Table 8. The 23 unpolarized and polarized SDMEs for p production from the deuteron in xb bins defined by the 
Hmits 0.0, 0.05, 0.08, and 0.35. The first uncertainties are statistical, the second systematic. 



Element 




{a;s)=0.064 


{a;B)=0.119 






0.410 ± 0.017 ± 0.032 


0.413 ± 0.012 ± 0.014 


0.411 ± 0.011 ± 0.010 




1 


0.233 ± 0.019 ± 0.018 


0.263 ± 0.013 ± 0.022 


0.223 ± 0.011 ± 0.012 


Im 


rU 


-0.226 ± 0.020 ± 0.039 


-0.234 ± 0.011 ± 0.019 


-0.218 ± 0.012 ± 0.018 


Re 


rfo 


0.162 ± 0.005 ± 0.021 


0.168 ± 0.005 ± 0.007 


0.148 ± 0.004 ± 0.011 


Im 


rfo 


-0.162 ± 0.003 ± 0.017 


-0.160 ± 0.004 ± 0.004 


-0.132 ± 0.005 ± 0.009 


Im 


rlo 


0.102 ± 0.038 ± 0.028 


0.106 ± 0.026 ± 0.006 


0.099 ± 0.025 ± 0.005 


Re 


rfo 


0.080 ± 0.017 ± 0.018 


0.080 ± 0.020 ± 0.005 


0.145 ± 0.024 ± 0.018 


Re 


r?^ 


0.044 ± 0.005 ± 0.009 


0.027 ± 0.005 ± 0.012 


0.026 ± 0.005 ± 0.005 


Re 


rlo 


-0.017 ± 0.013 ± 0.027 


-0.016 ± 0.010 ± 0.009 


-0.025 ± 0.009 ± 0.018 


Im 


rlo 


0.013 ± 0.008 ± 0.031 


0.024 ± 0.009 ± 0.019 


0.011 ± 0.010 ± 0.017 


'-00 




0.138 ± 0.014 ± 0.043 


0.075 ± 0.011 ± 0.014 


0.120 ± 0.010 ± 0.014 


'•oo 




-0.010 ± 0.023 ± 0.061 


-0.056 ± 0.022 ± 0.011 


-0.019 ± 0.024 ± 0.009 


Im 


rlo 


0.048 ± 0.026 ± 0.020 


0.010 ± 0.018 ± 0.005 


0.040 ± 0.017 ± 0.002 


'^00 




0.199 ± 0.043 ± 0.024 


0.074 ± 0.058 ± 0.005 


-0.059 ± 0.062 ± 0.016 


r!i 




0.010 ± 0.007 ± 0.021 


-0.006 ± 0.004 ± 0.005 


-0.030 ± 0.004 ± 0.015 


rl- 




0.018 ± 0.009 ± 0.011 


0.000 ± 0.005 ± 0.002 


0.020 ± 0.004 ± 0.007 


Im 




-0.010 ± 0.011 ± 0.002 


0.001 ± 0.005 ± 0.008 


-0.012 ± 0.005 ± 0.005 


Im 


ri-i 


-0.081 ± 0.033 ± 0.032 


-0.028 ± 0.037 ± 0.003 


-0.081 ± 0.036 ± 0.010 


rfi 




-0.028 ± 0.026 ± 0.007 


0.020 ± 0.029 ± 0.001 


0.049 ± 0.026 ± 0.004 




1 


-0.116 ± 0.042 ± 0.028 


-0.046 ± 0.035 ± 0.004 


-0.073 ± 0.034 ± 0.017 


r?i 


1 


0.004 ± 0.008 ± 0.020 


-0.003 ± 0.007 ± 0.012 


-0.007 ± 0.006 ± 0.002 


ri\ 




-0.023 ± 0.016 ± 0.020 


-0.007 ± 0.009 ± 0.011 


-0.010 ± 0.009 ± 0.010 


Im 


rl-i 


-0.001 ± 0.013 ± 0.019 


-0.032 ± 0.022 ± 0.004 


0.027 ± 0.019 ± 0.004 



Table 9. The values of the phase difference 5 between Tn and Too amplitudes calculated according to (|50p for the proton 
and deuteron in bins defined by the limits 0.5, 1.0, 1.4, 2.0, and 7.0 GeV^. The first uncertainties are statistical, the 
second systematic. 



target 


(Q^) = 0.82 GeV^ 


(Q^) = 1.19 GeV^ 


(Q^) = 1.66 GeV^ 


(Q^) = 3.06 GeV^ 


proton 
deuteron 


36.17 ± 12.33 ± 7.09 

33.18 ± 4.55 ± 9.88 


18.76 ± 6.99 ± 5.62 
22.55 ± 4.07 ± 2.83 


26.52 ± 3.92 ± 0.91 
30.37 ± 3.13 ± 0.54 


36.53 ± 2.79 ± 3.65 
36.58 ± 2.28 ± 3.99 



Table 10. Values in different kinematic bins of the variable ui = 1 — r-Qo + 2r?ii — 2r\-i — 2rJ_i, used for the test of 
NPE dominance, for proton and deuteron data. The first uncertainties are statistical, the second systematic. 



bin 


ui proton 


111 deuteron 


0.5 GeV^ < < 1 GeV^ 
1.0 GeV^ < 0^ < 1.4 GeV^ 
1.4 GeV^ <Q'' < 2.0 GeV^ 

2.0 GeV^ < < 7 GeV^ 


0.114 ± 0.053 ± 0.045 
0.148 ± 0.035 ± 0.044 
0.063 ± 0.037 ± 0.077 
0.178 ± 0.038 ± 0.040 


0.104 ± 0.035 ± 0.061 
0.069 ± 0.026 ± 0.048 
0.078 ± 0.028 ± 0.028 
0.169 ± 0.032 ± 0.024 


0.0 GeV^ < -t' < 0.04 GeV^ 
0.04 GoV^ < -t' < 0.10 GeV^ 
0.10 GeV^ < -t' < 0.20 GeV^ 
0.20 GeV^ < -t' < 0.40 GeV^ 


0.197 ± 0.043 ± 0.035 
0.090 ± 0.040 ± 0.041 
0.073 ± 0.041 ± 0.078 
0.125 ± 0.040 ± 0.107 


0.091 ± 0.029 ± 0.024 
0.082 ± 0.032 ± 0.039 
0.097 ± 0.033 ± 0.068 
0.077 ± 0.036 ± 0.095 


0.0 <XB < 0.05 
0.05 <xb < 0.08 
0.08 <XB < 0.35 


0.142 ± 0.099 ± 0.005 
0.123 ± 0.037 ± 0.029 
0.152 ± 0.031 ± 0.055 


0.180 ± 0.063 ± 0.075 
0.070 ± 0.027 ± 0.031 
0.149 ± 0.027 ± 0.039 
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Table 11. Results on iti, 112 and 113, calculated according to (|52)I54I) . shown together with average value or range in 
and W . Top section: HERMES results from proton and deuteron data, shown with statistical and systematic uncertainties 
separately. Bottom section: results from other experiments calculated from published SDMEs, with statistical and 
systematic uncertainties combined in quadrature without accounting for correlations between the SDMEs. 



Experiment 


Q^ GeV^ 


W, GeV 




U2 


U3 


HERMES V 


1.95 


4.8 


0.125 ± 0.021 ± 0.050 


-O.OUi 0.004 ± 0.012 


0.055 ± 0.045 ± 0.006 


HERMES d 


1.94 


4.8 


0.091 ± 0.016 ± 0.046 


-0.008± 0.003 ± 0.010 


-0.040± 0.035 ± 0.007 


ZEUS DIS [To] 


2.4 


90 


0.018 ± 0.066 


0.018 ± 0.011 




ZEUS BPC [9] 


0.41 


45 


0.058 ± 0.078 


-0.002± 0.016 




HI m 


2.5-3.5 


30-100 


0.065 ± 0.15 


-0.017± 0.034 




SLAC [55] 


0.9 


3.14 


0.85 ± 0.32 


-0.050± 0.072 




SLAC [56] 


0.9 


3.14 


1.174 ± 0.379 


0.039 ± 0.082 




DESY [32] 


1.05 


2 - 2.8 


0.73 ± 0.33 


-0.040± 0.064 





Table 12. Ratios of certain helicity-flip amplitudes to the square root of the sum of all amplitudes squared: roi 
for the transition 7^^ — > p'^, no for the transition 7^ — > Pt, and ti_i for the transition 7*^ — > p^. Top section: 
HERMES results from proton and deuteron data calculated according to (|59ll6ip . shown with statistical and systematic 
uncertainties separately. Bottom section: results from other experiments calculated according to (|62I64|I , with statistical 
and systematic uncertainties combined in quadrature without accounting for correlations between the SDMEs. 



HERMES 


TOl 


no 


Tl-l 


proton 
deuteron 


0.114 ± 0.007 ± 0.010 
0.122 ± 0.007 ± 0.006 


0.075 ± 0.030 ± 0.003 
0.090 ± 0.022 ± 0.011 


0.051 ± 0.029 ± 0.010 
0.007 ± 0.025 ± 0.015 


Experiment 


T"01 


TlO 


n-i 


ZEUS BPC|9] 
ZEUS DIS [10] 

HI m 

SLAC [55] 
SLAC [56] 
DESY [37] 


0.069 ± 0.027 
0.078 ± 0.016 
0.088 ± 0.036 
0.095 ± 0.165 
0.084 ± 0.177 
0.041 ± 0.247 


0.003 ± 0.029 
-O.OlOi 0.028 
0.019 ± 0.065 
0.030 ± 0.133 
0.412 ± 0.430 
0.335 ± 0.436 


0.048 ± 0.028 
0.013 ± 0.032 
0.035 ± 0.109 
0.112 ± 0.231 
0.042 ± 0.389 



Table 13. The longitudinal-to-transverse cross section ratios i?"*, iZ, and R^^^ for the proton and deuteron in 
bins defined by the limits 0.5, 1.0, 1.4, 2.0, and 7.0 GeV^. The total uncertainties are shown. 



Ratio 


Target 


(g^) = 0.82 GeV^ 


(Q2) = 1.19 GeV^ 


(Q2) = 1.66 GeV^ 


(Q^) = 3.06 GeV^ 




proton 
deuteron 


0.694 ± 0.187 
0.748 ± 0.179 


0.701 ± 0.063 
0.755 ± 0.063 


0.798 ± 0.080 
0.973 ± 0.098 


1.053 ± 0.074 
0.870 ± 0.061 


R 


proton 
deuteron 


0.649 ± 0.188 
0.677 ± 0.180 


0.671 ± 0.065 
0.583 ± 0.067 


0.794 ± 0.083 
0.958 ± 0.102 


1.068 ± 0.087 
0.922 ± 0.074 




proton 
deuteron 


0.755 ± 0.190 
0.809 ± 0.182 


0.783 ± 0.068 
0.798 ± 0.067 


0.840 ± 0.090 
1.041 ± 0.102 


1.225 ± 0.084 
0.994 ± 0.065 
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Table 14. The SDME results from the proton data, integrated over the entire HERMES kinematic range, presented 
in the notation of Ref. The first uncertainties are statistical and the second systematic. 





= 




0.412 


± 


0.010 


± 0.010 


Re - Uo+O) 


= 


V2{lm rto-rlo) 


-0.464 


± 


0.005 


± 0.028 


+ + 2eu++ 


= 




0.588 


± 


0.010 


± 0.010 




= 


7-1-1 - Im '-Li 


0.473 


± 


0.012 


± 0.029 


Reu°o% 


= 




-0.077 


± 


0.006 


± 0.006 


v°+ -11-° 4-2™°+ 


= 


2Re rO;^ 


0.062 


± 


0.008 


± 0.016 


Reii'i+ 




Re rJo - Im rfo 


-0.054 


± 


0.009 


± 0.027 


Re {n°~ - n+:^) 


= 


^(Im r-fo + Re rfo) 


-0.008 


± 


0.005 


± 0.015 


Re iu-+ + eUo+) 


= 


„04 
'1-1 


-0.011 


± 


0.005 


± 0.005 


Refill 




^11 


-0.025 


± 


0.007 


± 0.008 


Re + ^io+ ) 


= 




-0.023 


± 


0.004 


± 0.018 


Re 


= 


(Im rti-rLi)/\^ 


-0.005 


± 


0.005 


± 0.008 




= 


Tm 


0.011 


± 


0.019 


± 0.008 




= 


Re rJo + Im r-Jg 


0.009 


± 


0.009 


± 0.004 


Re Wg^T 




-(Im r?_i + rLi)/\^ 


-0.002 


± 


0.005 


± 0.001 


Re m1+ 




rj_i + ImT-J„i 


0.018 


± 


0.012 


± 0.011 


Im («°+ - u^°) 




Im rio + Re rfo 


0.264 ± 0.030 ± 0.023 


Im 




r!o/y2 


0.025 ± 0.035 ± 0.007 


Im («•]+ - 




-2Im r?o 


0.034 ± 0.030 ± 0.008 


Im (^°; - ^++0) 




y2(Im rIo - Re rfo) 


0.054 ± 0.035 ± 0.008 


Im 




— Im ri_i 


0.024 ± 0.018 ± 0.001 


Im + 






0.51 ± 0.034 ± 0.001 






(rti + Im rLi)/\^ 


-0.012 ± 0.039 ± 0.007 






(rtl-ImrLl)/^/2 


0.038 ± 0.039 ± 0.001 



SOME 


„04 
' 00 


^04 


„04 
'^1-1 






rlo 


rl 

'^1-1 


rlo 


r2 

'^1-1 


r-fi 


'■oo 


rfo 


^5 

'^l-l 


'^lO 


ri-i 


rlo 


'^1-1 


rlo 


r^ 

'^1-1 


rfi 


'^00 


„8 
'^10 


„8 

'1-1 


'^00 


1.00 












































0.16 


1.00 










































-0.05 


-0.01 


1.00 






































r\. 


-0.04 


-0.01 


0.62 


1.00 




































rlo 


-0.02 


0.12 


-0.18 


-0.48 


1.00 


































Re rJo 


0.00 


-0.41 


-0.08 


-0.03 


0.02 


1.00 
































rLi 


-0.65 


-0.22 


-0.05 


0.00 


0.00 


0.02 


1.00 






























Im rlo 


0.02 


0.33 


-0.03 


-0.03 


0.07 


-0.17 


-0.10 


1.00 




























Im rf_i 


0.54 


0.23 


0.00 


-0.03 


0.01 


-0.11 


-0.32 


0.07 


1.00 


























r?i 


-0.17 


0.15 


0.02 


-0.04 


0.09 


-0.03 


0.10 


0.07 


-0.07 


1.00 
























'-oo 


0.37 


0.01 


0.06 


0.09 


-0.26 


0.02 


-0.27 


0.00 


0.20 


-0.50 


1.00 






















Re rfo 


-0.10 


0.08 


-0.13 


-0.12 


0.06 


-0.18 


-0.03 


0.08 


0.04 


-0.01 


-0.12 


1.00 




















rti 


-0.05 


-0.15 


0.10 


0.10 


-0.09 


0.01 


0.04 


0.02 


-0.08 


-0.43 


0.08 


-0.04 
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